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XVII. BAND SECHSTES (SCHLUSS-) HEFT 1950 


Die Biegeschwingungen eines Stabes unter pulsierender Achsialkraft 


bei beliebigen Randbedingungen. 
Von K. Kloiter. 


1. Aufgabenstellung. Die Biegeschwingungen, die ein Stab unter der Wirkung pulsierender 
Achsialkrafte ausfiihrt, sind schon mehrfach behandelt worden!. Das Problem wurde danach in 
_allgemeinere Zusammenhinge gestellt? und schlieBlich zu einer Theorie der Stabilitit elastischer 
Bewegungen itberhaupt ausgebaut®. Daneben wurden an Einzelproblemen noch untersucht die 
Biegeschwingungen eines Kreisringes unter pulsierendem radialem Druck und die Kippschwin- 
gungen eines [-Tragers bei pulsierenden Endmomenten?. 

Fiir den Stab unter pulsierender Langskraft (das Problem, von dem alle diese Untersuchungen 
ihren Ausgang nahmen), ist eine explizite Lésung bisher allerdings nur in dem einen Sonderfall 
bekannt, in dem der Stab an seinen beiden Enden gelenkig (momentenfrei) gelagert ist. Die sonst 

_ noch méglichen Randbedingungen schienen das Problem bisher so zu komplizieren, daB eine Lé- 
sung nicht gefunden werden konnte. Es laBt sich indessen zeigen, und das soll in dieser Mitteilung 
_ geschehen, daB auch diese Faille einer unmittelbaren Behandlung zuganglich sind. 
Um den Vergleich der neu in die Untersuchung einbezogenen Falle mit dem schon bekannten 
_ Fall des beiderseits gelenkig gelagerten Stabes méglichst tibersichtlich zu machen, wird im Ab- 
schnitt 2 eine nochmalige Diskussion dieses Falles (mit einer z.T. andersartigen Darstellung der 
~ schon bekannten Tatsachen) vorausgeschickt. Es folgt im Abschn. 3 die Integration der Be- 
wegungsgleichung unter den neuen Randbedingungen; daran schlieBt sich im Abschn. 4 eine Dis- 
kussion der so gewonnenen Lésungen hinsichtlich ihrer Stabilitat an. Der letzte Abschn. 5 faBbt 
die neuen Ergebnisse noch einmal zusammen. 
2. Der beiderseitig gelenkig gelagerte Stab. Die Differentialgleichung der Querbewegung 
y (x, t) eines unter einer pulsierenden Achsialkraft P= P)-- S cos Qt stehenden Stabes mit der 
(konstanten) Biegesteifigkeit EJ und der Langendichte « lautet, wie man aus der Beziehung 
Mo = EJ (M Biegemoment, 0 Kriimmungshalbmesser) leicht herleitet°, 


EJy!¥ + (P+ Scos Qt) y+ wy =0. (1) 
Sie ist eine partielle, lineare Differentialgleichung, die hinsichilich der Langenkoordinate x von 
vierter, hinsichtlich der Zeit t von zweiter Ordnung ist. Ihre Koeffizienten sind nicht alle kon- 
stant;im zweiten Glied tritt als Koeffizient die Funktion P(t), eine Funktion der Zeit, auf. Der 
Produktansatz 


\y 


y (x, t) = Y(x)- Ti), (2) 

der bei konstanten Koeffizienten stets (und deshalb auch noch bei unveranderlicher Achsial- 

- kraft P) zu einer Trennung der Veranderlichen verhilft (so daB man fiir die Funktion Y (x) eine 

gewohnliche Differentialgleichung vierter Ordnung, fiir die Funktion T(t) eine gewohnliche Diffe- 

rentialgleichung zweiter Ordnung erhalt) bewirkt diese Trennung bei der Differentialgleichung (1) 
nicht mehr. 

Wenn der Stab an beiden Enden gelenkig gelagert ist, wenn also die Funktion y (x, t) den Rand- 
bedingungen 


y(O,)=0, ¥"(0,2)=9, yLj=% yy (L)=0 (3) 
zu allen Zeiten t geniigen soll, so 1aBt sich mit Hilfe des besonderen Ansatzes 
y (x, t) = A sin (max/l)- T(t) (2a) 


1 E. Mettler, Mitt. Forsch.-Anst. Gutehoffnungshiitte 8 (1940), 5. 1. — K. Klotter, Forsche. Ing.-Wesen 12 
(1941), S. 209. — E. Mettler, Forschungshefte Geb. Stahlbau H. 4 (1941), 5.1. 
2 EF. Meitler, Ing.-Arch. 13 (1942), S. 97. 
8 E. Mettler, Ing.-Arch. 16 (1947), S. 135 sowie 17 (1949), 5. 418. 
4 S. Woinowsky-Krieger, Ing.-Arch. 13 (1942), 5. 90 u. 197. , 
5 K. Klotter, a.a.O. Gl. (22) mit geringer Abanderung der Bezeichnungen. 
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wegen des in allen drei Gliedern der Differentialgleichung (1) in gleicher Weise erscheinenden und 


deshalb herausziehbaren Faktors sin (mz x/l) allerdings die Abhangigkeit von der Ortskoordinate | 


« beseitigen und eine gewohnliche Differentialgleichung fiir die Zeitfunktion T(t) gewinnen. Diese 
Differentialgleichung Jautet, nachdem man noch die Transformation Qt = T vorgenommen hat, 


eT 1 /EJ mint — P, m2? min? es fs 
daa Fal nw I a 12 > ann or seat ll ( ) 

Sie ist eine Mathieusche Differentialgleichung vom Typ 
T’+ Ti+ y cos t)=— 0 


mit tT als unabhangig Veranderlicher. Unter Benutzung der Abkiirzungen 


(5) 


, EJ mx‘ (Eigenfrequenzquadrat der Biegeschwingungen m-ten Grades des un- 
Te belasteten Stabes, P,= 0), 
Pa jaye (Knicklast des Stabes mit m Halbwellen), (6) 
ir Pid ake AES 
ie like : . Prim | 


erhalt die Differentialgleichung (4) die Fassung 


2 


(3) 
eee Le Pres) eo Cog t/t. —=0r 


Aus ihr entnimmt man die Parameter der Differentialgleichung (5) fiir die Funktion T(t), 
w? or, 
= pl (5a) 
Hinsichtlich der Funktion T(t) interessiert nicht so sehr ihr Verlaufim einzelnen als vielmehr 
die Frage, ob sie beschrankt bleibt (stabiler Fall) oder nicht (instabiler Fall). Hinen Uberblick 
iiber jene Wertepaare A, y der Parameter, die zu stabilen oder zu instabilen Lésungen T(t) fiihren, 
y gibt die Stabilitatskarte 
der Mathieuschen Differential- 
gleichung, die sog. Ince-Strutt- 
sche Karte (auch Struttsche 
s Mire Karte schlechthin genannt) 
VEN WS LA (Abb.1). Beschrankt man die 
N Betrachtung weiterhin auf den 
vornehmlich interessierenden 
SSNs Fall m= 1 (eine Halbwelle) 
So : (wir lassen im folgenden den 
Zeiger 1 bei w, P,, p und s 
einfach weg), so kann man aus 
der Stabilitatskarte (Abb. 1) 
& SS die Diagramme der Abb. 2 
< \ gewinnen, in denen jeweils 
die (mit P, dimensionslos ge- 
machte) Schwankungsampli- 
tude S als Ordinate, der (eben- 
so behandelte) Mittelwert P, 
als Abszisse aufgetragen ist, 
wahrend der Quotient @/Q 
von Diagramm zu Diagramm 
andere Werte hat. Die Dia- 
gramme zeigen, da, wenn eine von Null verschiedene Schwankungsamplitude S vorhanden ist, 
1. stabile Bewegungszustande auch dann noch vorhanden sein kénnen, wenn der Mittelwert 
P, den Eulerschen Wert der Knicklast, Pg, iiberschritten hat, 
2. dafs auch bei Mittelwerten P,, die unterhalb des Wertes Pg liegen, ja sogar beim Mittelwert 
Py = 0, instabile Bewegungszustinde auftreten kénnen, Fri . 
3. dafs sogar im Bereich negativer Werte P), d. h. wenn die Langskraft P, eine Zugkraft ist, 
Solche Instabilitaten vorkommen kénnen. 


(4a) 


m 
—— oO Sm > 


| | C, 
| rs 


Abb. 1. Ince-Struttsche Karte fiir die Mathieusche Differentialgl. (5). Schraffierte 
Gebiete: stabil; helle Gebiete: instabil. 


1 K. Klotter, a.a.O. Gl. (25). 
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‘ i Ce Weise findet man die Tatsachen in den bisherigen Darstellungen angefiihrt. Statt 
er es a Serie von Diagrammen nach Abb. 2, in denen die GréBen s und p als Koordi- 
naten dienen*, verwendet E. Mettler ein gemeinsames Diagramm, bei dem als Ordinate der Quo- 


Abb, 2. s,p-Stabilitatskarten der Funktionen r() (t) fiir verschiedene Frequenzverhiltnisse. 
(a) w?/Q? = 1/8, (b) w2/Q2 = 1/4, (c) w9/Q? = 3/4, (d) w9/Q2=1, (e) w2 QI= 4. 
tient 2/q@,, (in unserer Bezeichnungsweise), als Abszisse der Quotient S/(m? Pr, -— Py) auf- 
getragen ist. 


Wir wollen jetzt zunachst noch eine Vereinfachung vornehmen. Erstens bedenken wir, daB 
die Schwankungsamplituden S in der Regel klein gegen die Knicklast P, sein werden, so daf wir 


la 
oY 


N 

t= N 
ae NX! 

@ IS 
as SSX 
14 x oN 

i Si 
12 == = NS SN 


S 0 LY | INS: 4 a 
‘ VAIL SW 
FE NNNSURIE 


iS 
FY SNNNNEES 
ADDO MAA AWMTZANS MYON 
FES IRRRAA RRA 


ANN 
SX 
S 


Nt 


-0§ -0% -03 -02 -07 U0 Of 2 0607 +08 09 10 


Abb. 3. y,A-Stabilitatskarte bei Anwesenheit von Dimpfung. Lage der Crenzkurven fiir verschiedene 
Dampfungswerte 0; 107-4; 107%; 107°; Si dOses LOS, 
mit Werten s <1 rechnen diirfen. Fiir solche kleinen Werte s sind wegen der durch (5a) aus- 
gedriickten Proportionalitat zwischen s und y die Instabilitatsbereiche schmal, wie ein Blick auf 
Abb. 1 lehrt. Im Hinblick auf die (hier zwar explizit nicht beriicksichtigte, aber doch stets vor- 
handene) Dampfung, durch die die Spitzen der Instabilitatsbereiche von der A-Achse aus hoch- 


1 Siehe FuBnote 1 von S. 363. 
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gezogen werden, und zwar um so mehr, je héher die Ordnung (angegeben durch den Parameter A) 
des Gebietes wird (Abb. 3)!, ist itberdies praktisch allein der erste Instabilitatsbereich (dessen Spitze 


im Punkte 7 = 


tigung finden. Beschrankt man die Betrachtungen auf die Schwingungsform mit einer Halbwelle 


= liegt) von Bedeutung; allenfalls mag der zweite noch gelegentlich Beriicksich- 


(m = 1) und auf den ersten Instabilitatsbereich, also auf Werte Ain der Nahe von = , so gibt die 


erste Gleichung (5a) den Zusammenhang zwischen dem Mittelwert Py der Langskraft und ihrer 
Schwankungsfrequenz 2 an, der zu instabilen Lésungen fihrt. Mit der Abkiirzung 


Q 
ane (7a) 
lautet er 
2 
P= P,(1—") (7b) 


und wird (fiir den Bereich der Druckkrafte, also positiver Werte P,) durch die Parabel der Abb. 4 
verdeutlicht. Fiir Wertepaare Q, P,, die Punkten ,,in der Nahe‘: dieser Kurve entsprechen, hat 


p 
1A, PE2 
Pe 
S2 
7-0; 
| ee “Fy 
a 7 2 ett 
Abb. 4. Die mittlere Druckkraft P, in Ab- 1-& 
hangigkeit vom Frequenzverhialtnis 2/a,. 
Tnstabilitatskurve fird = 1/4 und m= 1 beim o h a 3 . 6 6 4 a 
gelenkig gelagerten Druckstab; Abb. 5. Weitere Instabilitatskurven fiir den pulsierend belasteten 
Funktion r*)(2). Druckstab; Funktionen 7 *) (1) fiir m=1 und m=2. 


man also unbeschrankt anwachsende Funktionen T(t) zu erwarten. Besonders gefahrlich sind 

solche Lingskrafte, die mit einer Frequenz 2 = 2 @, pulsieren; denn dann fithren schon Krafte 

mit dem Mittelwert Null zur Instabilitat der Biegeschwingungen. 
In Abb. 5 sind auBer der in Abb. 4 gezeigten (fiir 


4A, 1 
d i Fe und m= 1 geltenden) Kurve noch diejenigen, 
i S analog zu verstehenden, Kurven eingetragen, die sich auf 
das erste Instabilitatsgebiet (2 He =) aber die Lésung mit 
zwei Halbwellen m= 2 (ebenfalls ausgezogen) und ferner 
g Mdiejenigen, die sich auf das zweite Instabilitatsgebiet 
al 4 q To (A= 1 mit m—1 und m= 2, gestrichelt) beziehen. 


Die Angabe ,,in der Nahe‘+ kann dann prdazisiert 
werden, wenn man genauere Aussagen iiber den Betrag 
der kleinen, die Schwankungsamplitude S messenden 
GréBe s macht. Fiir s <1 findet man z. B. zu beiden Seiten der Parabel (7b) in Abb.4 die um 


Abb. 6. Instabilititsbereiche ,,in der Nahe“ 
der Kurve aus Abb, 4. 


1 
= S nach unten und oben verschobenen Parabeln (Abb. 6) 


2S 
Py= Pe| seek. (Ze) 


die die Grenzen des Instabilitatsbereichs angeben. Die Gleichungen (7c) folgen aus den Glei- 
chungen der die Grenzkurven c,;, und s,;. in Abb. 1 ersetzenden Tangenten, 


a= LNA (7d) 


durch Einfithren der GréBen Py, Pe und 7 gemaB den Abkiixzungen (6) und (7a). Ganz ent- 
sprechend sehen die Instabilitatsgebiete um die tibrigen in Abb. 5 stark ausgezogenen (zu A = re 


1 Vel. FuBnote 1 von S. 363 sowie G. Kotowski. Z. angew. Math. Mech. 23 (1943), S. 226. 
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gehérenden) Parabeln aus, wahrend die Instabilitatsbereiche um die gestrichelten (zu A= 1 ge- 
hérenden) Parabeln erheblich schmaler bleiben, sofern sie nicht vermége der Dimpfung iiber- 
haupt verschwunden sind. 


3. Integration der Bewegungsgleichung unter anderen Randbedingungen. Die Differential- 
gleichung der Querschwingungen ist in jedem Fall die Gleichung (1). Liegen aber Randbedin- 
gungen vor, die von (3) verschieden sind, so gibt es keinen Ansatz vom Typ (2) oder (2a) mehr, 
der erlaubte, die Abhingigkeit von der Ortskoordinate x dadurch zu beseitigen, daB ein gemein- 
samer Faktor aus allen Gliedern der partiellen Differentialgleichung (1) herausgezogen wird, so 
da eine gewéhnliche Differentialgleichung fiir die Zeitfunktion iibrig bliebe. Wegen dieses Um- 
standes scheiterten die bisherigen Bemithungen zur Integration der Differentialgleichung (1), 
wenn andere Bedingungen als (3) vorlagen. Eine kleine Abwandlung der Uberlegungen eréffnet 
hier aber trotzdem einen Weg zum Ziel. 

Im Ansatz (2a) ist die Kreisfunktion sin(mz«/l) die Eigenfunktion des an beiden Enden ge- 
lenkig gelagerten querschwingenden Stabes; sie erfiillt die Randhedingungen (3). Die Eigenfunk- 
tionen anderer Lagerungsfalle sind, wie immer auch die Randbedingungen im einzelnen lauten 
moégen, Aggregate aus den vier Kreis- und Hyperbel-Funktionen 

COS Hm Xs SIN Hm Xs Co} %m x» Sin 4m x (8) 
mit Koeffizienten, die die Kreis- und Hyperbelfunktionen des Argumentes %,1 enthalten; die 
GréBe x, geniigt dabei der jeweiligen Frequenzengleichung 

(och, t) == 0. 3 
Als Beispiele betrachten wir 
a) den an beiden Enden fest eingespannten Stab, fiir den also die Randbedingungen 
y(O.)=0, y(0)=0, y(hy=0, y= (92) 
bestehen, und 

b) den an einem Ende gelenkig gelagerten, am andern Ende fest eingespannten Stab mit den 

Randbedingungen 
¥ (054) = 0; yA { 0; 2) == il bpt) ==0,, yi (Lt) == 0. (9b) 

c) Zum Vergleich schreiben wir die Bedingungen (3) des schon behandelten, an beiden Enden 

gelenkig gelagerten Stabes als Gleichungen (9c) noch einmal an: 
w102) == 0, y= (0, t)' == 0, a (Ld) == 0"; Val t= Or, (9c) 

Die Eigenfunktionen Y*(x) der freien Querschwingungen der Stabe lauten in den drei Fallen 

der Reihe nach 
Yn, (x) = (cos %ml — Co] xml) (cos %m x — Co} %m x) 


+ (sin %ml+- Sin xml) (sin %m x — Sin %m x), (10a) 
YS (x) = Sin xml sin %m X — Sin Xml Sin %m Xx ; (10b) 
5 il (Cg es TOW ee ae (10c) 
als Frequenzengleichungen treten auf 
1 — cos %ml Co} xml = 0, (lla) 
tS nl 29 nl == 0, (11b) 
sine) 0 (11c) 


Damit eine Lésung y(x, t) der Differentialgleichung (1) anderen als den Randbedingungen (3) 
angepakt werden kann, muB sie aus Gliedern bestehen, die die Funktionen (8) midglicherweise 
alle enthalten. An Stelle des Ansatzes (2a) verwenden wir deshalb jetzt den vervollstandigten 
Ansatz 
y (x, t) = A 08 Hm %° TO (t) + B sin xm x* T)(t)+- C Cof xm %* T(t) + D Sit xm x: T(t) (12) 
mit zundchst noch unbestimmten Zeitfunktionen T(t). 

Gehen wir mit dem Ansatz (12) in die Differentialgleichung (1) ein und ordnen sogleich nach 
den Funktionen (8), so erhalten wir 


A008 Hm %{{EJ 4, — 22, (Py + S cos Q1)] TY + TH} + 
+ B sin %m x {{EJ xt, — x2, (Po S cos Q2)] T?) + wT} + (13) 
+ CGof xm x {LE J wn + xm (Po+ S cos 2t)] TS) + pw TO) 
+ D Sinem x{{E J xm + xm (Po + S cos 21)] T) + uw TH) = 0. 
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Die Gleichung (13) ist demnach erfiillbar, wenn die beiden Zeitfunktionen T(t) und T‘?) (t) der 
Differentialgleichung 


[EJ 24, — x2, (Py + S cos 2t)] T™ + w T® = 0, (14a) 
die beiden Funktionen T(t) und T(t) aber der Differentialgleichung 
[EJ xt, + x2, (Py + S cos Q)|T™ + wT = 0 (14b) 


geniigen. (Die oberen Zeiger (k) und (h) sollen dabei auf die Kreisfunktionen und Hyperbelfunk- 
tionen hindeuten.) Die beiden Differentialgleichungen (14a) und (14b) sind Mathieusche Diffe- 
rentialgleichungen. Beide lassen sich unter Benutzung der wegen der abweichenden Rand- 
bedingungen sinngemaB zu verstehenden Abkiirzungen @, Pr, m, p und s schreiben als 


2 
TH” + *1(1 — Pm) — Sp 008 7] TH = 0 (15a) 
und 
2 
TA Ut Pe) city CoB cr] EO a0 (15b) 


Eine partikulare Lésung y,, (x,t) der Differentialgleichung (1) lautet dann gemaB dem Ansatz (12) 
V¥m(X, t) = (Ap CO8 %m x + B,, sin %m%) Te 
+ (Cn Cof %m% + Dm Git xm %) Tm’ (t)- 


Wegen der Linearitat ist auch jede Summe solcher Lésungen wieder eine Lésung, so daB die 
allgemeine Lésung lautet 


(16a) 


y (xt) = 3) ym (x8). (16b) 


n= 


e 


4, Diskussion der neuen Lésung. Die Koeffizienten A,, bis D,, in der Lésung (16b) miBten 
nun so bestimmt werden, daB die Randbedingungen, beispielsweise die Bedingungen (9a) oder 
(9b), erfiillt werden. Diese Aufgabe erforderte die Auflésung eines Systems unendlich vieler al- 
gebraischer Gleichungen der Art, wie sie hier beispielsweise gemaB den Bedingungen (9a) an- 
geschrieben werden: 


Sed eo) Ste) Gy TEN 
SR BET OO = = DT Gye 
Dy (Ag C08 %m_1 + By sin %, 2) re (t) = — S'(C,, Cof xm 1+ Dm Sin xml) Tee (t), - (17) 


Dy (— %m Am Sin %m 1+ Xm Bm COS %m!) r©) (t) 


Alle Summen sind dabei iiber m von 1 bis oo zu erstrecken. Das Gleichungssystem (17) stellt 
vier Gleichungen dar zwischen cot Konstanten, Gleichungen, die zu allen oo! Zeitpunkten t be- 
stehen miissen. 

Zur angenaherten Lisung des Problems innerhalb eines vorgegebenen Zeitabschnitts kénnte 
man sich darauf beschranken, in endlich vielen, etwa n, Zeitpunkten t,, t....t, die Gleichungen 
(17) zu erfiillen. Man erhielte so 4n Gleichungen, die zur Ermittelung von jeweils n der vier Grup- 
pen von Konstanten A,, bis D,, dienen kénnten. 

Da uns aber, wie wir schon mehrfach betonten, nicht so sehr der zeitliche Ablauf der Bewe- 
gung im einzelnen interessiert als vielmehr deren Stabilitat, so verzichten wir auf die Ermittelung 
der Integrationskonstanten und schlieBen folgendermafen weiter: Die Lésung y(x, t) der Diffe- 
rentialgleichung (1) hat die Gestalt (16b). Damit sie beschrankt bleibt, miissen alle Summanden 
¥m(%,t) und somit alle Funktionen TG) und 7 (t) beschrankt bleiben. Diese Funktionen 
sind aber Lésungen der beiden Mathieuschen Differentialgleichungen (15a) und (15b). Die erste, 
die Differentialgleichung fiir die Funktionenklasse T (t) ist formal identisch mit der Differen- 
tialgleichung (4a), die wir im Abschn. 2 schon eingehend erértert haben. Die Differentialgleichung 
(15a) unterscheidet sich von der Differentialgleichung (4a) nur dadurch, daB fiir sie die Zeichen 
Om PE, m und damit p,, und s,, eine allgemeinere Bedeutung haben. Die GréBen w,, und Piss 
bezeichnen nicht Eigenfrequenz und Knicklast des Lagerungsfalles (3) sondern des jeweils vor- 
liegenden Lagerungsfalles. Unter Beriicksichtigung dieser Wandlung in der Bedeutung der Zei- 
chen gelten die fritheren Aussagen und somit auch die Diagramme der Abb. 2, 4 und 5 fiir alle 
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Lagerungsfalle, wenigstens hinsichtlich der Partikularlésung m —1. Man muB allerdings beach- 
ten, daB zur Erfillung von Randbedingungen, die von (3) verschieden sind, neben der Partikular- 
lésung m =1 auch die hiéheren herangezogen werden miissen. 

Die zweite Funktionenklasse T (t) geniigt der Differentialgleichung (15b). Interpretiert 
man sienach genau denselben Gesichtspunkten wie die Lésungen der Differentialgleichung (15a), 


WN 


LLL 


Abb. 7. s, p-Stabilitatskarten der Funktionen 7") (t); (a )bis (e) wie in Abb. 2. 


so kommt man statt zu den Diagrammen der Abb. 2 zu denen der Abb. 7 (sie sind einfach an der 
Ordinatenachse gespiegelt) und statt zur Kurve der Abb. 4 zu der der Abb. 8. Die Instabilitats- 

A, punkte (P,,.Q), dieentweder durch die Lésungen 
T®) oder die Lésungen T“*) bedingt werden, 
liegen demnach (im Bereich der Druckkrafte, 
P, > 0) in der Nahe der Kurven der Abb. 9. 


a a, 20; D2 2O> 
ae a ; 4 Abb. 9. Instabilitatskurven beim pulsierend belasteten Druck- 
Abb. 8. Instabilitatskurve fiir 2 = 1/4 und m = Acres pulsie stab bei Beriicksichtigung beliebiger Randbedingungen; 
rend belasteten Druckstab; Funktion 7! Y(t.) Funktionen rik) (¢) und Te), 
m m 


5. Zusammenfassung. Wahrend fiir die Instabilitat des an beiden Enden gelenkig gelagerten 


Stabes die Instabilitatsgebiete der der Differentialgleichung (14a) geniigenden Funktionen To 
allein maBgebend sind, miissen in allen anderen Lagerungsfallen auch die Instabilitatsgebiete der 
der Differentialgleichung (14b) geniigenden Funktionen T beriicksichtigt werden. Die In- 
stabilitatspunkte P,, Q liegen (fitr Kleine Werte s) in der Nachbarschaft der Parabeln der Abb. 9. 
Die Nachbarschaft kann bei Kenntnis der Schwankungsamplitude S fiir die zum ersten Instabili- 
tatsbereich (4 —1/,) gehérenden (ausgezogenen) Parabeln so abgegrenzt werden, wie Abb. 6 


andeutet. 
Diese Aussagen gelten, ganz gleichgiiltig wie die von (3) verschiedenen Randbedingungen im 


einzelnen aussehen, also nicht nur fiir die als Beispiele herangezogenen Bedingungen (9a) und 
(9b), sondern fiir alle méglichen. 


(Eingegangen am 5. Juli 1949.) 
Anschrift des Verfassers: Prof. Dr.-Ing. K. Klotter, Karlsruhe-Riippurr, Hegaustr. 4. 
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Flache kegelige Korper in Uberschallstr6mung. 
Von H. Behrbohm und K. Oswatitsch. 


1. Begriff und Theorie der kegeligen Strémung wurde von A. Busemann' gegeben. Flache 
kegelige Stérkérper in homogener Uberschallparallelstrémung lassen eine Linearisierung der be- 
schreibenden Differentialgleichungen zu. Im Bereich 1,5 << M, <5 der Machzahlen der An- 
strémung wird dadurch das Wesentliche bereits wiedergegeben. In der vorliegenden Arbeit 
werden die friiheren Ergebnisse auf neuem Wege hergeleitet. Busemanns Transformation auf 
die Laplacesche Gleichung und Heranziehung der konformen Abbildungen im Strémungsbereich 
innerhalb des Machkegels mit Spitze im Kegeligkeitszentrum wird vermieden, statt dessen findet 
Pucketts Methode der Belegung mit quellartigen Singularitaten? Anwendung. Da bei kegeligen 
Feldern die Zustinde auf konzentrischen Strahlen konstant sind, lassen sich die flachenhaften 
Integrale der Puckettschen Methode durch Ausfihren einer Integration auf Integrale tiber eine 
Variable reduzieren. Das gegeniiber der bisherigen Behandlungsweise Neue der vorliegenden 
Note ist die Tatsache, daB sich diese einfachen Integrale durch geeignete partielle Integration 
auf die Form Poissonscher Integrale bringen lassen. Dadurch wird das Umstrémungsproblem 
flacher kegeliger Kérper — bei jeder Lage dieses Kérpers relativ zum Machkegel durch seine 
Spitze — mittels der singularen Integralgleichung der Wirbelschicht beschrieben. Diese Integral- 


gleichung spielt sowohl in der Prandilschen Theorie der tragenden Linie als auch in der Theorie * 


diinner Tragfliigelprofile eine zentrale Rolle. Thre Lésung ist unter weiten Voraussetzungen be- 
kannt?. 

An Literatur seien zusatzlich genannt die Arbeiten von Stewart*, Sauer> und Germain‘, die 
alle den von Busemann eingeschlagenen Weg beschreiten, und die Arbeiten von Brown’ und 
Ferri®, die sich der Puckettschen Methode unter Hinzunahme von Dipolen bedienen. 


2. Im Raum der rechtwinkligen kartesischen Koordinaten x’, y’, 2’ befinde sich ein flacher 
kegeliger Stérkérper St’ mit Spitze im Ursprung 0, der seine stationare Parallelanstrémung in 
x'-Richtung (Machzahl M. >1, Anstrémgeschwindigkeit u’,) nur wenig stért. Die resul- 
tierende Strémung werde als isentrop und drehungsfrei angenommen, und uw’, v’, w’ seien ihre 
Geschwindigkeitskomponenten in x’-, y’-, 2’-Richtung. Es werden nun folgende Gré®en ein- 
gefiihrt : 


ae y=¥ (M1, z= VME I, | 
u=(W—uL)VME—1, =", w= w' | 


Die Isentropieforderung gibt die Dichte 9 als Funktion des Druckes p. Die Tatsache, daB §’ 
flach ist und also nur wenig stért, erlaubt, auf approximativem Wege die Kontinuitatsgleichung 
zu linearisieren, indem Glieder zweiter und héherer Ordnung in den Stérgeschwindigkeiten u, v, w 
vernachlassigt werden. Die Drehungsfreiheit gestattet die Einfiihrung eines Stérpotentials 
p(%, ¥, 2) mit p,— u, py = v, p, = w. Fiir dieses y geht mit Hilfe der genannten Linearisierung 
die gasdynamische Gleichung in die dreidimensionale Schwingungsgleichung 


Px — Pyy — Paz = 0 (2) 
itber. Zu dieser Differentialgleichung treten noch die Randbedingungen u= v= w= 0 im 
Anstrémgebiet und die Randbedingung am Storkorper, die zum Ausdruck bringt, daB die Normal- 
komponente der Geschwindigkeit auf die Kérperoberflache verschwinden muB. 


Die Oberflache von §t, des Bildes von St’ zufolge (1), werde durch F(x, y, z) = 0 beschrieben. 
Das Wesentliche dieser Abbildung besteht u. a. darin, daB der Machsche Winkel der Anstrémung 


(1) 


x eis Busemann, Schriften der Deutschen Akademie der Luftfahrtforschung Bd. 7B (1943), H. 3, 


* A. E. Puckett, J. Aeronautical Sciences 13 (1946), S. 475. 
° H. Séhngen, Math. Z. 45 (1939), S. 245. 
* H. J. Stewart, Quart. appl. Mathematics, USA, Okt. 1946, S. 246. 


| R. Sauer, Berechnung der aerodynamischen Beiwerte der Leitwerke nicht rotierender Geschosse bei 
Uberschallgeschwindigkeit, unveréffentlichtes Manuskript. 

6 P. Germain, Office National d’Etudes et de Recherches Aéronautiques, Nr. 34, 1949. 

PC, 15e Brown, N. A.C. A. Technical Note No. 1183 (1946). ig 

8 A. Ferri, Elements of aerodynamics of supersonic flows. New York 1949, S, 350—383. 
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45° wird. DaB §& flach sei, werde nun etwas praziser dahingehend formuliert, daB F, und (i 
Klein seien gegen F,. Unter diesen Annahmen kann aus der Normalkomponentenforderung 


F,(u+u.)+Fy,o+F,w=0 fir E'(e./45;2) == 0 


mit Hilfe der Forderung gegen ux kleiner Stérgeschwindigkeiten u,v, w auf die Naherungs- 
gleichung 


Pepe oy =O) fur F(x... 8) =. 0 


geschlossen werden, welche Gleichung ihrerseits dann wieder statt auf der Kérperkontur 
FP(«, y, z) = 0 auch auf z= 0 gefordert werden kann. Im Rahmen der fiir die vorgenommene 
Linearisierung der Strémungsgleichungen zu treffenden Vernachlassigungen kann daher die 
Umstrémungsforderung von § in einfacher Weise als geforderte Vorgabe einer w,-Verteilung 
Wy = w(z= 0) in demjenigen Gebiet G der (x, y)-Ebene angesehen werden, auf das sich ® in 
2-Richtung projiziert. G ist infolge der Kegeligkeitsforderung ein durch zwei von 0 ausgehende 
Halbstrahlen begrenzter Winkelraum (oder eine Summe von solchen; doch bietet ein solcher 
allgemeinerer Fall fiir das Ergebnis keine anderen Schwierigkeiten als solche der Darstellung, 
weshalb er hier nicht behandelt werden soll). 

Die Eigenschaft der Kegeligkeit hat zur Folge, da8 der Strémungszustand, also u, v, w, nur 
von y/x und z/x abhangen darf; u,v, w miissen also homogene Funktionen nullter Ordnung, 
d.h. p(x, y, z) mu homogene Funktion erster Ordnung von x, y, 2 sein. 


3. Der EinfluBkegel des Ursprungs 0 (Kreiskegel mit halbem Offnungswinkel von 45° und 
mit positiver x-Achse als Achse) heiBe auch der Machsche Kegel It der Strémung. Das Innen- 
gebiet von t werde als das elliptische Gebiet E der Strémung bezeichnet. Es hat namlich mit 
den Gleichungen von elliptischem Typ die Eigenschaft gemeinsam, in allen seinen Teilen von der 
Anderung einer kegeligen Lésung betroffen zu werden. Denn jede solche Lisung reicht in die 
unmittelbare Umgebung des Ursprungs, deren EinfluBgebiet eben durch It gegeben ist. Das 
Gebiet H auBerhalb Ii, abgesehen von dem von 0 aus stromaufwarts symmetrisch zu ‘it ge- 
legenen Kegel, in welchem kegelige Lésungen nicht auftreten kénnen, werde als hyperbolisch 
bezeichnet. E, bzw. H, sei der Durchschnitt von E bzw. H mit der Ebene z= 0. E, besteht aus 
den inneren Punkten des ersten und achten Oktanten, H, aus den inneren Punkten des zweiten 
und dritten Oktanten einerseits, des sechsten und siebenten Oktanten der Ebene z= 0 an- 
dererseits. 

Hadamards Grundlisung von (2), die dreiparametrige Funktion 


1 4 
ped fan fir x— x > (y— %)? +(2—*)?, 
P(x: ¥, 2) = | |@—*?—-G—yF—-e—aP Vy ya) ) 


0 sonst, 


ist der mathematische Ausdruck fiir das Geschwindigkeitspotential einer quellartigen Singularitat 
im Punkte x, ¥o, 2) bei linearisiertem Uberschall. Indem man jeden Punkt x, 0, 0 der positiven 
x-Achse mit einer solchen Quelle des Gewichtes —x, belegt, erhalt man gema dem Super- 
positionsprinzip bei linearen Differentialgleichungen die Lésung 


eV ee sae 
‘ | eS ha ogg a eee rae Pa 
g,(P) = He Cat a pea eee 
| 0 fir P auBerhalb E 


von (2). Sie ist homogen erster Ordnung in x, y, 2, entspricht also einer kegeligen Strémung. 
Pist hier und im folgenden der Aufpunkt mit den Koordinaten x, y, z. Durch analoge Belegung 
der positiven y-Achse wird die Lésung 


y+ Ve _—- 
Yo4 Yo — ¥(5 are sin —— \—j#—¥—# fir Pin E, 
\#?—(y—y)? —# 2 ye—2?, 
olP)=| os 
a ‘ Yo 4 Yo = ny fiir Pin H 
, L420.) 2 
y— \x?—28 


erhalten. 
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Bei —1 < m <1 geht (2) mittels der ,,Lorentztransformation“ 


, B= %,9(%, ¥, 2) = (4%, Y> 2) (3) 


in die Differentialgleichung yy. — 455 — 4a = 9, d.h. in sich ttber. Demgema8 wird aus 7, (P) 
sofort eine Lésung y,(P; m) von (2) erhalten, wenn x, y, z statt der dortigen x, y, 2 geschrieben 
und diese dann gema® (3) durch x, y, z ausgedriickt werden. g,(P;m) entspricht einer der 
x-Achsenbelegung analogen Quellbelegung der in E, gelegenen Halbgeraden y= mx, x > 0, 
|m| <1. Da St bei (3) in sich tibergeht und ein in E oder H gelegener Punkt auch nach der 
Transformation in E oder H bleibt, sind Belegungen von in H, gelegenen Halbgeraden auf diese 
Weise nicht zu erreichen. Vielmehr hat man dazu von Y(P) auszugehen und den analogen Uber- 
gang gem4B (3) zu machen wie bei y,(P). Man kommt so zur Lésung ~,(P; n) von (2), wenn zur 
Unterscheidung mim elliptischen, nim hyperbolischen Bereich fiir den Transformationsparameter 
aus (3) geschrieben wird. g,(P; 1) entspricht einer der (positiven) y-Achsenbelegung analogen 
Quellbelegung der in H, gelegenen Halbgeraden x= ny, y > 0, |n| <1. 

Analog hatte man auch die negative y-Achse belegen und aus ihr durch hyperbolische Drehung 
die Belegungen von in y < 0 gelegenen Halbgeraden von H, gewinnen kiénnen. Diese in jedem 
konkreten Fall leicht zu treffenden Erganzungen sollen zur Vereinfachung der Darstellung jedoch 
nicht mit aufgefiithrt werden, zumal sich im folgenden das Interesse nur auf Belegungen von E, 
und den in y > 0 gelegenen Teil von H, konzentrieren wird. : 


4. Durch abermalige lineare Superposition werden aus y,(P;m) und 9,(P; 7) allgemeinere 
kegelige Lésungen ©, und @®, von (2) gebildet, indem ,(P; m) bzw. y,(P;n) mit den von den 
Richtungskoeffizienten m bzw. n des jeweils belegten Halbstrahls abhangigen Gewichten e(m) 
bzw. h(n) multipliziert und danach iiber den ganzen Bereich —1 < m <1 bzw. —l<n<l 


nach m bzw. n integriert werden: 


1 


D,(P)= [ 9,(Psm) e(m)dms — @,(P) =f y(Psm) h(n) dn. (4) 


Von nun an interessieren nur noch die durch die Quellbelegungen hervorgerufenen Stér- 
geschwindigkeiten, d.h. die partiellen Ableitungen der @;(i—= 1,2) nach x, y,z. Da diese 
homogene Funktionen nullter Ordnung von x, y, z sind, ist es bequem, die homogenen Koordi- 
naten 7 = y/x, = 2/x einzufiihren. Zundchst sei nur die z-Komponente w(7,¢) der E,-Be- 
legung betrachtet. Man ermittelt fiir ¢ 40 die Darstellung 


1 


——————— (1 — m?) e(m) dm 3 ; x 
= Soe 2 4 2 2 
w (7, ¢C) Pet. n t it need 2) fir Pin E, d.h. fir 7? +2 <1, 


ai! 


wahrend w fiir P auBerhalb E identisch verschwindet. Hieraus folgert man durch partielle 
Integration [unter Annahme der Existenz von lim (1 — m2) e(m)] leicht das Puckettsche Er- 
gebnis eee 


Wo(7) = w(n, 0) ie w(y, C) = a(1 — 7) e(n) fir alle —l<7y< 1. (5) 
+ +0 
Es gibt einen sehr einfachen Zusammenhang zwischen der Gewichtshelegung e(m) der elliptisch 
gelegenen Halbgeraden und der von dieser Belegung auf dieser Halbgeraden erzeugten z-Kompo- 
nente w)(m) der Stérgeschwindigkeit, wobei, was wohl beachtet werden mége, der Grenziiber- 
gang ¢—> --0 aus dem oberen Halbraum £ > 0 heraus durchzufiihren ist. Im Falle ¢— — 0, 
d.h. von ¢ < 0 her, hatte sich nach dem dritten Gleichheitszeichen in (5) ein Minuszeichen er- 
geben, was fiir spater wichtig ist. 


Indem man nun e(m) in (4) durch w,(m) ausdriickt und die Funktion 


SS ais Soi hoa Pere al ees ye 
f(m) = f(y, 6; m) = ——] l fiir 72 2 
iim) ee 
einfiihrt, kann man die aus der E,-Belegung herrithrenden Stérgeschwindigkeiten u(7, 6), v(n, 6), 
w(7,C) bei Beschrankung auf den oderen Halbraum in der Gestalt 
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1 
Ce) .¢3;m 
u(y, ¢)= —— i ge ) Wy(m) dm, 
—i 
1 

1 Of (yn,c3;m 
| — ) yey(m) dm, fir Pin E } (6) 

—i 

i il 
OS eo | ane (ets 4 (rm dm | 

/ ST) ge 


ah 


schreiben. Fiir €—> +0 geht (6) in 


(4) = w (0) = Him (5 2) = — [OEY om) dn 
= (7) 
V (1) = v (7, 0) = lim v (7, 6) = 


o++0 


1 
a em 2) Wy (m) dm 
=il 

fir Pin E, d.h. fir —1 < 7 <1 itber, wo itber die singulare Stelle m= 7 des Integranden 
hinweg hier und im folgenden mittels Cauchyscher Hauptwertbildung zu integrieren ist. Um 
die Stérgeschwindigkeiten u,v des durch (4) gegebenen Stérpotentials O,(P) auch im unteren 
Halbraum richtig zu bekommen, hat man der obigen Bemerkung zufolge ne rechten Seiten in 
(6) und (7) die entgegengesetzten Vorzeichen zu geben, da die erzeugten u, v ja gerade Funk- 
tionen von € sind. 

AuBerhalb E, d.h. auf und auBerhalb 9, verschwinden die von ®,(P) herrithrenden u, v, w 
identisch. 

Ganz das Entsprechende gilt nun auch fiir die aus ®,(P) d. h. aus der H,-Belegung herkommen- 
den Stérgeschwindigkeiten. Zunachst gilt analog zu (5) wieder das Puckettsche Ergebnis 


Wy (n) ue win 6) = + a(1 — n?) h(n) fir alle —-l<n<l. 


Fiir ein hyperbolisch und im EinfluBgebiet der Belegung gelegenes P gelten dann bei Beschran- 
kung auf den oberen Halbraum und unter Benutzung der Abkiirzung 


_ xy —|z|Vy?+2— 2 
om yee 


die Formeln 


2Q 
)d + — x? wy (2) 
al ea Ne Ne a xz We tea = % : ate 
2 —_—_—_ — ——_— —— 
wy(n)dn yz (y? + 2? —2 x*) + a (P—y?) VY +H — 2? w(Q) 
v(x, y, 2) = (i — ns? (72 + 2) 10 


= 


2 Piet )\\ Q) 
w(x, ¥, j= tree =) jeer 


die fir z— + 0 in 


x/¥ 
NW, » )dn Wo (x/¥) _ 
U(x, ¥) = (1 — n2)3?2 on: ’ 
=e 
x/ 
PEGS dn, % wy (x/y) 


g(x, ¥) = — (1 — nn)? Tees 


iibergehen. 
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Zur einfachen Darstellung der in einem elliptisch gelegenen P durch die H,-Belegung bewirkten 
Stérgeschwindigkeiten fithrt man die Funktion 


g(n) = 2(y, C3 nn) = + are sin — ips - for 9-6" |n| << Il 


yl a VQ —nn)? — (1 — n?) 2? 


ein. Mit ihrer Hilfe findet man bei Beschrankung auf den oberen Halbraum die Formeln 


ail Da 
uOnd= Z f EGE mlm) an, 


TU 


] 


» (nC) = ale (ng (7 65) y(n) dn, fiir Pin E (8) 


On 


bzw. 


1 
Oso 
uy (n)= + [lee ) wy (n) dn, 


Tt On 
vie & fiir Pin E (9) 
il 0 BO ) 
% (n)=— 4 i RUSE FAC is! 


—1) 


Im Gegensatz zu (7) sind die Integrale in (9) gewéhnliche Integrale. Im unteren Halbraum haben 
die Formeln fiir u und v die entgegengesetzten Vorzeichen. 


5. Entscheidend fiir den hier vorzulegenden neuen Aufbau der Theorie linearisierter kegeliger 
Strémungen ist nun eine Umformung der Formeln (6) bzw. (7) fiir die Stérgeschwindigkeiten u, v 
aus einer E,-Belegung. Man ziehe aus den Integranden der ersten beiden Formeln (6) den Faktor 
1/(1 — m?)$/? heraus und integriere alsdann partiell. Hierzu bilde man aus der Geschwindig- 
keitshelegung w)(m) die Integrale 


a Wo d d 
Fy(m)= [ Gene und F,(m) = fe ee ~ (10) 


fiir |m| <1 und fihre die beiden Funktionen 


K,(m)= me a (eee eu F,(m)) = aut am __ _F,(m) 2 | 
Fees (11) 
__ %(m) F,(m) | 


— m2 V1—7: mae 


a [m| <lein. Mit diesen Bezeichnungen gehen dann aus (6) durch partielle Integration die 
ormeln 


1 
mw u(y, C) = V1— 72 — £2 m— 1 OP ee w,(1) see 
une) = V9 on jo 


1 
Reyes {_ 72 m—n 7 2 Le ( Mo(1) wo(— 1) 


fiir 7,¢ aus E und € £0 hervor. Der Grenziibergang € > -L0 liefert hieraus 
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1 
aE Se K,(m) dm 1+) 1/1 —7 
a(n) = 17 | m—1 (V3 = ep (1) & = 5 1). 
1 


J fiir —1 <4 <1. \(12) 


jl—7 (Sas (2 


"2 
/ m—y 1—y7 


gt 0,,(7) = w, (1) pee Wo ( 1) 


Dabei sind die nach Anwendung der Hopitalschen Regel auf die ausintegrierten Bestandteile 
erhaltenen Ausdriicke w (1) und w)(— 1) als die Grenzwerte von wy(m) bei m—> 1 — 0 baw. 
bei m—> —1 +-0 (also vom Innern von E, her) einzusetzen. Gleichung (12) gibt die in der 
Einleitung bereits erwahnte Poissonsche Integraldarstellung der von einer E,-Belegung her- 
riihrenden Stérgeschwindigkeiten in der (x, y)-Ebene. 

6. Mit ®,(P) und ®,(P) aus (4) ist auch O(P) = O,(P) + @,(P) Lésung von (2). Sie er- 
fillt offenbar die Randbedingungen u = v = w= 0 im Anstrémgebiet. Damit @(P) das Stor- 
potential der Uberschallumstrémung des kegeligen Stérkérpers § sei, ist nach den Entwick- 
lungen des Abschnitts 2 noch die Bedingung 


FF, 
W(x, ¥) == —U,, r fir x, yin G 


zu erfiillen. § wird fiir zs > 0 durch Vorgabe von w, auf dem oberen Ufer von G, fiir z < 0 durch 
Vorgabe von w, auf dem unteren Ufer von G gegeben. Fiir jeden der beiden Halbraume kénnen 
dann mit dem Formelapparat aus 4 bzw. 5 die Stérgeschwindigkeiten wu, v, w fiir z< 0 durch 
Wo, U, v, w fiir z < 0 durch w, ausgedriickt werden. Danach sind diese beiden Strémungsbereiche 
auBerhalb G den Randbedingungen gema8 zusammenzufiigen. Dabei sind zwei Falle zu unter- 
scheiden: 

I. wy und w, sind auf der ganzen Ebene z= 0 bekannt. Dies ist einerseits der Fall, wenn 
beide Rander von § in H liegen, andererseits wenn § bei beliebiger Lage seiner Rander symme- 
trisch in bezug auf z= 0 ist. Stets ist dann w= wy = 0 auBerhalb CG. 

II. wy und w, sind in einem Teil von z= 0 zunachst unbekannt. Dies ist dann der Fall, 
wenn ein asymmetrischer Kérper das Gebiet E, nicht ganz erfiillt. Auf Gist dann w, und wy, vor- 
geschrieben, auf E, auBerhalb G muB w= w, sein. 

Nur dieser Fall [list noch zu behandeln. Er la8t sich durch Lorentztransformationen auf zwei 
spezielle Falle zuriickfihren: 

1. Ein Korperrand ist die x-Achse, der andere liegt in H, und ist gegeben durch x= ny, y, 
z= 0 fir y > 0 und |n,| < 1. 

2. Beide Kérperrander liegen symmetrisch zur x-Achse in E, und sind durch die Halbgeraden 
y= + mx, z= 0 fir x >0 und 0 < m, < 1 gegeben. 

Im Unterfall 1. ist w, und w, fiir G gegeben und fiir E, auBerhalb G, d.h. fir —1 << m < 0, 
so zu bestimmen, da$ dort u,(m) = u,(m), v)(m) = v,(m), w)(m) = wy(m) gilt. Man zeigt mittels 
der Gleichung du,/dy = 0v,/@x fiir die Wirbelfreiheit leicht, daB die v 9-Gleichheit keine neue 
Bedingung darstellt, sondern stets gleichzeitig mit der u)-Gleichheit erfiillt ist. Fiir einen Auf- 
punkt P(7) mit —1 < 7 < Oist nach (9) und (12) unter Beachtung von w,(— 1) = w)(— 1) = 0 

1 0 
1+ ——= { K,(m)dm ee a K,(m) dm 
Et ntl) +yi—a [ See" + yim | BS 
, U 


1 


miig(q) = | 20%”) ran) an — J/ 


ny 


aus dem oberen Halbraum her, nebst einer analogen Gleichung aus dem unteren Halbraum her, 
in der nur rechts jedes Glied das entgegengesetzte Zeichen tragt. In den rechten Seiten beider 
Gleichungen sind jeweils nur die Hauptwertintegrale unbekannt. Aus der Gleichsetzung der 
linken Seiten erhalt man daher unter Beachtung von K,(m) = K,(m) = K,(m) fir —1 << m< 0 
die singulare Integralgleichung 


f | 
K,(m) d ‘ 
An = Lyn) fir —1l<n<0 
—] 


mit bekannter rechter Seite L,(7). Sie hat die Lésung (Fufn. 3) 
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0 
C 1 1 | L,()V—tt+)) g, (13) 


K = SS 
(mm) 22 /— m(m + 1) Te" J—m (m+ 1)_ m—t 


=) 


fiir —1 < _m < 0, mit zunichst noch unbekannter Konstanten C. Hieraus kann mittels (10) 
und (11) auch w,(m) fiir —1 < m < 0 angegeben und durch dieses Zuriickgehen auf w,(m) auch C 
bestimmt werden. Damit ist dann das in G bereits gegebene w,(m) und w,(m) zu einer auch in 
ganz E, definierten Funktion erganzt, also schlieBlich in der ganzen Ebene z= 0 bekannt. 
Im Unterfall 2. ist (1) = ™(—1) = wo(1) = w)(— 1) = 0 und daher fur ein P(y) mit 
Many <1 nachs(12) 
—m, 1 


cjlath- ) nem sia) Ketan wlempl==ees fir) hme eee (ad Re ea 
nig) = ie | Sede pyre | SO yim | Se 


m— 1 
—J] m, mL 


mit analoger Gleichung fiir 7u)(7), in der dann rechts nur lauter Minuszeichen stehen. Aus der 
ip-Gleichheit folgt daher 


—m, 1 an LUD 
, [4 [FESO aml —— [EOVLBOD a tie |r| < ful <] 
aay m, 


mit bekannter rechter Seite. Infolge der Kérpersymmetrie in bezug auf y — 0, die hier lediglich 
der Einfachheit der Darstellung und ihres vorzugsweisen Auftretens in der Praxis wegen voraus- 


gesetzt werden soll, ist K,(—m) = — K,(m), K,(—m) = — K,(m), so daB man mit der Abkiirzung 
K,(m) + K,(m) = 2 K,(m) auf die Integralgleichung 


rx, (Va) a 


t— 7" 


m2 

— i AE fiir, m2 << 
m? 0 

mit bekannter rechter Seite gefiihrt wird. Dies ist wieder eine Integralgleichung vom selben 
Typ wie im Unterfall 1., aus der K,(m) und damit w,(m) fiir |m,| << m < 1 gefunden werden 
und damit dann das gesamte Strémungsfeld berechnet werden kann. 

7. Mit Hilfe der oben dargelegten Entwicklungen kann das réumliche Geschwindigkeitsfeld 
um irgend einen flachen, kegeligen Stérkérper in homogener Uberschallparallelstrémung im 
Rahmen der in Abschnitt 2 geschilderten Linearisierung der Strémungsgleichungen berechnet 
werden. Derselbe Kérper § wird danach bei zwei verschiedenen Anstellwinkeln durch eine um 
den konstanten Betrag der Anstellwinkeldifferenz geanderte w - und w,-Verteilung gegeben. 
Mit anderen Worten: Ist die Druckverteilung des nicht angestellten Kérpers bekannt, so wird 
seine Druckverteilung bei irgendeiner kleinen Anstellung dadurch aus jener gewonnen, da man 
ihr die Druckverteilung der entsprechend angestellten, unendlich diinnen ebenen Platte additiv 
iiberlagert. 

Wegen der hierdurch aufgedeckten groBen Rolle, die die diinne tragende Dreiecksplatte fiir 
die Berechnung angestellter kegeliger Kérper spielt, und weil ihre Theorie als reife Frucht leicht 
aus den bisherigen allgemeinen Darlegungen entnommen werden kann, seien zum AbschluB die 
wichtigsten Angaben iiber diesen an und fiir sich bestens bekannten Gegenstand erlaubt. 

Fir das elliptisch (und symmetrisch zur x-Achse) gelegene Dreieck der Abb. 1 ist ww (m) = 
wo(™m) = w)(m) = A= const fir —m,<m<m,, daher K,(m)= K,(m)= K(m)=—0 fir 
—m, << m< m,, so daB die rechte Seite von (14) verschwindet. Man hat die Lésung 


Kin) e filr,| it, <n <0 


V1 — mi? mn? — m? 


mit K,(—m) = —K,(m). Hieraus ergibt sich nach (10) und (11) 


wy(m) = A+ C ee i dt | 


2 a 1 ——— 
i iG ili 238 ||P Se 
J 8yi—#je—m 


und also wegen w,)(1)= 0 


m? A+ CE'(m,) = 0, 
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wo E"(m,) das zum. Modul = m+ gehorige vollstandige Legendresche Integral zweiter Gattung 
bezeichnet. Damit ist C bestimmt. Nach elementaren Rechnungen ergibt sich aus (12) 
ee 1 a n 
u,(n) = 3 —, (Nn) = 2 — = 
o( ”) | E’ (m,) Vm? — — ne U 0 (77) oe tie (m,) | m2? — 7? 
fir —m, <7 <m,. Damit sind samtliche drei Stérgeschwindigkeiten an der Dreiecksfliche 


ermittelt. Die oberen Vorzeichen gelten fiir das obere, die unteren fiir das untere Ufer der 
tragenden Flache. 


Fiir das gemischt ge- 
legene Dreieckder Abb. 2 
ist 


W)(m) = wo (m) = Wy (m) 
== A —=.const 
fur Ocerm =] 
und 
Wo (n) aS Wo (7) = w,(n) 
==, = Gn 
fit tee mt =< | 


bei —l cn,z<l. 
Daherist K, (m)—= K,(m) 
== (fir 0 < m < 1 und 


2 Abb. 1. Das elliptisch gelegene, zurx-Achse Abb,2. Das gemischt gelegene Dreieck mit zur 
fiir das L, der rechten symmetrische, tragende Dreieck. ‘Anstrémung paralleler Randkante, 
Seite von (13)ergibt sich 

A 


If = 
(7) hay eae 
Damit wird nach (13) 


& +haresin 7) fire’ ly 0 


Il 10.77 


Kn) = Bose = eae esl ace et 

a V1 —n? V1 —m? Jl—n, /—m+ /1—m 

wenn die Konstante C jener Formel durch Zuriickgehen auf w)(m) aus der Forderung der Endlich 
keit des Massenflusses geeignet bestimmt wird. Man erhalt so aus (10) und (11) 


| mie == Nl ein< WU. 


2A Ve A a le aot ri ee 
uo (m) = — Tae | — a parcsin TOE m)m) fir —l<m<0 
und daraus nach (12) 
eT ne pone Lh CRSA i 
OOS ara ea ar ae funk Oey) <a 
und durch Vermittlung von K,(m) schlieBlich noch 
ee ee ee ene Pin = fir O<y <1. 
v(y)=+ = eee F ta are sin |! Ui 


Alle anderen Lagen tragender Dreiecke werden durch Lorentztransformationen auf die beiden 
hier behandelten zurickgefiihrt. 


(Eingegangen am 6. Marz 1950.) 


Anschrift der Verfasser: Dr. H. Behrbohm, Waldkirchi. Breisgau, Merklinstr. 38 
Dozent Dr. K. Oswatitsch, Stockholm, Flugtechn. Inst. d. Techn. Hochschule. 
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Beitrag zu einer Wirbeltheorie des kompressiblen Mediums. 
Von J. Dorr. 


1. Aufgabenstellung. Wir wollen das von einem bewegten Schallstrahler oder von einem im 
Fluge schwingenden Fliigel erzeugte Geschwindigkeits- und Druckfeld des umgebenden gas- 
férmigen oder fliissigen Mediums untersuchen. Dabei wollen wir annehmen, da8B sich das Problem 
in zwei Raumdimensionen beschreiben lasse. In einem kartesischen Koordinatensystem (x, y) 
erfolge die Bewegung des Schallstrahlers in Richtung der negativen x-Achse. In Richtung der 
y-Achse gemessen sei der Schallstrahler so flach, das man ihn in der (x, y)-Ebene durch ein auf 
der x-Achse liegendes Geradenstiick ersetzen kann. Die Schwingungen des Schallstrahlers seien 
beschrieben durch zeitabhangige Deformationen und Verschiebungen dieses Geradenstiickes in 
y-Richtung. Diese seien so klein, daB man bei der Darstellung des Geschwindigkeitsfeldes des 
umgebenden Mediums Gréfen von héherer als erster Ordnung vernachlassigen kann. Dann darf 
man das Problem mit der Differentialgleichung der Schallausbreitung 


1 ao 
behandeln. Dabei ist ¢ die Schallgeschwindigkeit und 
v, = grad, » (1b) 


die Geschwindigkeit des Mediums in Richtung des Vektors 3. Die durch die Differentialgleichung (1) 
gestellte Aufgabe mu noch erginzt werden durch die Randbedingung, daB das Medium den 
schwingenden Schallstrahler immer liickenlos umschlieBen muS8. Wenn die Marschgeschwin- 
digkeit v des Schallstrahlers kleiner als die Schallgeschwindigkeit cist, dann mu zur Gewahrung 
einer eindeutigen Lisung noch die Kuttasche AbfluBbedingung beachtet werden. Wir betrachten 
im folgenden nur den Fall v < c. 

Fiir ¢—> oo geht (1a) in die klassische Potentialgleichung itber, die fiir zwei Dimensionen mit 
funktionentheoretischen Methoden behandelt werden kann. Bei instationdren Problemen fiihrt 
aber die AbfluBbedingung zu einer wesentlichen Erschwerung der gestellten Aufgabe. Trotzdem 
konnte fiir den Sonderfall c— oo mit funktionentheoretischen Methoden eine exakte Lésung 
gefunden werden und zwar von Ellenberger-Schmieden! und von Theodorsen?. 

Kin anderer Weg zur Lésung bietet sich fiir c—> 00 mit einer auf Birnbaum zuriickgehenden 
Methode; bei der die Lésung durch eine Uberlagerung von gebundenen und freien Wirbeln ge- 
wonnen wird. Dieses Verfahren fiihrt auf eine recht komplizierte Intelgralgleichung, deren exakte 
Lésung Schwarz? und Séhngen4 gelang unter Heranziehung des von Betz fiir den stationadren 
Fall gefundenen lésenden Kerns. 

Eine Ubertragung dieser Methoden auf den Fall des kompressiblen Mediums ist nicht ohne 
weiteres méglich. Die funktionentheoretische Methode versagt, weil sich die Differential- 
gleichung (la) fiir 0? g/t? +0 nicht funktionentheoretisch deuten last. Die Wirbelmethode 
versagte, weil fiir das kompressible Medium kein Potentialausdruck bekannt war, mit dem sich 
ein gebundener Wirbel mit seinem am Entstehungsort liegen bleibenden freien Wirbel (Anfahr- 
wirbel) beschreiben last. Prandtl léste an sich alle Schwierigkeiten durch Einfiihrung des sog. 
Beschleunigungspotentials, bei dem sich an Stelle der Elementarlésung des Wirbels die des 
Dipols verwenden la8t. Die auBerordentliche Fruchtbarkeit der Wirbelmethode fiir das inkom- 
pressible Medium sollte aber ein Hinweis sein, die Ubertragbarkeit der Wirbelmethode auf das 


kompressible Medium in aller Konsequenz zu iiberpriifen. Das Folgende soll ein solecher Versuch 
sein. 


2. Der gebundene Wirbel konstanter Auftriebskraft mit seinem freien Anfahrwirbel im kom- 
pressiblen Medium. Wir fiihren unsere Betrachtungen zunichst durch in bezug auf ein karte- 
sisches Koordinatensystem (x’, y,z), das bezogen auf das ungestérte Medium ruht. Die ein- 


+ Ellenberger und C. Schmieden, Z. angew. Math. Mech. 16 (1936), S. 193. 
2 Theodorsen, NACA-Bericht Nr. 4.96. 

3 L. Schwarz, Luftfahrt-Forsch. 17 (1940), S. 379. 

4 Séhngen, Luftfahrt-Forsch. 17 (1940), S. 401. 
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fachste Elementarliésung der Differentialgleichung (1a) ist das Potential einer ruhenden Punkt- 
quelle. 


pp =— f(et—r) mit r= V(x’ — x1)? + (y — 4)? + (2— 2). 


Darin bedeutet f eine beliebige beschrankte Funktion. Das Quellzentrum liegt bei (x,, ¥,, 2,). 
Kine Linienbelegung solcher Quellpunkte ist ein Quellfaden. Wir denken uns im Zeitpunkt 
t= 1, einen beidseitig unendlich langen Quellfaden konstanter Stirke 1/2 auf der z-Achse 
aufgelegt. Das zugehérige Potential hingt nicht von z ab. Mit 


Gea a(t —1,)? — x tye 


findet man dafiir den Ausdruck 


oo z an | PO (ip caeeee a OR, AD) 
ented | pa yt IEP 
2 J epee jeep 


gy (x’,y,t—t, =0 fir x2? +y?> ce? (t—1t,). 
Daraus findet man durch Differentiation das Potential eines Dipolfadens. 
og’ = ye (t —t,) 
Ge pute ala Whee sed ea cai 
Pp (x',y,t—t)=0 fir x2 +y2> ce? (t—1t,)?. 
Wir gehen jetzt auf ein Koordinatensystem itber, das sich in bezug auf das ungestoérte Medium 
mit der konstanten Geschwindigkeit v<c in negativer Richtung der x’-Achse bewegt. Es sei 


also 
x= x’ +v(t—t,), 
Po (%, ¥, t— ty) = Pp (x — v (t— ty), y, t— h). 

Es ist demnach pp das Potential eines ebenen Dipols, der im Zeitpunkt tt, spontan an der 
Stelle (x= 0, y= 0) entsteht und dann in bezug auf das Koordinatensystem (x, y) mit der 
Geschwindigkeit v nach rechts abschwimmt. 

Wir denken uns jetzt im Zeitpunkt tt anfangend laufend solche Dipole im Koordinaten- 
ursprung entstehen und nach rechts abschwimmen. Das Potential dieser Dipolbelegung ist 


Pp (x', ¥,t—t,)= 


Pw (X,Y, t— Tt) =v [Yo (%, yt — ty) dty. 


Da gp in einem gewissen AuBenbereich identisch Null ist, findet man mit der Substitution 


E=v (t == t,) 
und der Abkirzung 
2 , 
Beate aig ga AS EO) 2) 
den Ausdruck 
v (t—rt) 


yé _ dé 
Pw (%, Ys t—T) = ———— & 
(eer + yA |/P— Sle + 9 


a 


Das Integral ]a8t sich berechnen. Es ist, wie man durch Differentiauion nach 7 = v (t — 7) leicht 


bestatigen kann, 


yr @— ye — Sx — 0 OP +99 


= 3a) 
sical y+ afx—v(t—7)] ; ( 
m=0 fir (c?— v*) (¢— 71)? +2 0% (t—7)— — xy? <0. (3b) 
Fiir das inkompressible Medium erhalt man daraus fiir t 27 
yv(t—7) ys ian Te) 4) 
lim @,, = are tg ei oa eee ae are tg 7 are tg 7 ( 


c—>x 


27 
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Das ist der bekannte Potentialausdruck fiir einen am Orte (x= 0, y = 0) gebundenen Wirbel 
und seinem zugehérigen mit der Strémungsgeschwindigkeit v fortschwimmenden Anfahrwirbel. 
Der Ausdruck (3a) ist also das Analogon dazu fiir das kompressible Medium. Eine Aufspaltung 
von (3a) in zwei Summanden ist nicht wie in (4) derart méglich, da der eine Summand nur den 
gebundenen, der andere nur den freien Wirbel darstellt. Im Kompressiblen ist also das Wirbel- 
paar eine untrennbare Einheit. 

Wir wollen noch die Potentialfunktion 9,, fiir sehr kleine Entfernungen von den beiden Wirbel- 
zentren angeben. Dazu entwickeln wir den Wurzelausdruck in (3a) zuerst fiir sehr kleine Werte 


von x und y: 


[eX — vy yu 
yvr—r|/ Cc eee Tee: } 
Pw = are tg ES Sar 
Sere et eee ees fiir jety< {v (t—T) 


“| aa by eae | (Ew) (fsa) 


Der erste Summand stellt einen gebundenen Wirbel dar, der, wie es der Fall sein mu8B, nach der 
Prandtlschen Regel verzerrt ist. Der zweite Summand ist der Potentialanteil am Orte des 
gebundenen Wirbels, der vom freien Wirbel herriihrt. Fir kleine Entfernungen vom Zentrum 
des freien Wirbels findet man, wenn man alle unwesentlichen Terme fortlaft, 


Ge= — are tg <—"E—) +--- fiir Viz— v@t—)P +y?<v(t—tT). (5b) 


Das ist der Potentialausdruck fiir einen nicht verzerrten freien Wirbel. 
Fiir die vorgesehene Randwertaufgabe miissen wir noch die auf der x-Achse vorhandene 
Geschwindigkeitskomponente w, in y-Richtung kennen, die von dem Wirbelpaar erzeugt wird: 


re | @—2— Slee oF 


le ae = ol OES aM 2 (6) 


Ww, (x) ity fiir Wc? (t—1)?—[x— v(t P< 0. 


Wir sehen, dafi der ,,Abwind“* zwei Pole erster Ordnung aufweist, von denen einer am Orte des 
gebundenen, der andere am Orte des freien Wirbels liegt. 

Wir wollen noch die am Orte des gebundenen Wirbels wirkende Auftriebskraft berechnen. 
Dazu benétigen wir die Druckfunktion //, die den Druckunterschied zwischen dem 6rtlichen 
Druck und dem Normaldruck im Unendlichen angibt. Es gilt 


II (x, 9, )=—e(3, +34] Qu (7a) 


Darin ist @ die Dichte des ungestérten Mediums. Wir setzen zur Abkiirzung t— 7 —1t’ und 
erhalten 


: (7b) 


IT=0 fir c?t’*—(x—vt’')P—y?< 0. 
Die Druckfunktion besitzt, wenn man von dem Grenzkreis des Bereiches 
ce? t? — (x— vt’)? —y > 0 


absieht, nur am Orte des gebundenen Wirbels (x = 0, y =0) eine Unstetigkeit. Es kann des- 
halb nur an dieser Stelle eine Auftriebskraft wirken. Wir bezeichnen diese Auftriebskraft mit 
A, und berechnen sie in iiblicher Weise wie folgt: 

+a 


A, = apm, f [— IT (x, y, t) + II (x, —y, t)] dx. 


Sesh Ope ae) 
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Darin ist a eine beliebige, aber feste positiv reelle Zahl. Mit x — ul\y| findet man 
+ co 

/ 2 

i — Z = 


A,=20v = du=22700, (8) 


ice Va 
a 


See Che 


A,=0 fir t<r. 


Nachdem das Wirbelpaar entstanden ist, wirkt also am Orte des gebundenen Wirbels eine Auf- 
triebskraft, die unabhangig von der Lebensdauer des Wirbelpaares ist. 


3. Der gebundene Wirbel variabler Auftriebskraft mit seiner freien Wirbelschleppe im kom- 
pressiblen Medium. Wir nehmen jetzt an, es entstehen laufend als Funktion der Zeit t Wirbel- 
paare Pw», wie sie durch (3) beschrieben sind. Die Belegungsfunktion sei g(t). Der Vorgang 


beginne im Zeitpunkt t= — T. Als Potential dieser Wirbelbelegung findet man, da ~, in einem 
gewissen AuBenbereich identisch Null ist, 


pet 

@,, (x, y,t) = if 2(T) Dw (x, ¥,t — 7) dr, (9) 
“6 a 

G(x, ¥,t) =O tir as rt Sa 0. 


Die GréBe a ist durch (2) gegeben. Auf Grund von (8) kénnen wir sofort den am Orte des ge- 
bundenen Wirbels (x0, y= 0) wirkenden Auftrieb A des Wirbelsystems (9) angeben: 


A(t) =2z0v (ae dt. (10) 


Auch die Geschwindigkeitskomponente K in y-Richtung auf der x-Achse kénnen wir mit (6) 
sofort angeben: 
gee 


K(«,t)= f g(x) w,(x,t—1) dr. (11) 
7 
Hier ist nach (2) 
z vce|x| 
sh aa a a moe (12) 


Der Auftrieb A (t) ist eine beliebige Funktion der Zeit, da auch g(t) eine beliebige (integrierbare) 
Funktion ist. Das Wirbelsystem, welches im Punkte (x0, y—0) diesen Auftrieb liefert, wird 
durch (9) beschrieben. Fiir die vorgesehene Randwertaufgabe miissen wir insbesondere die 
durch (11) dargestellte senkrechte Geschwindigkeitskomponente auf der x-Achse kennen. Setzt 
man in (11) den Ausdruck (6) ein, so erhalt man mit t—t=6& 


t+T felt Uist Spee . 
K(n)= f e¢—2- 55) 1—-() (13) 


K(%,t)=0 fir v(t+T)< a). 

Dieser Ausdruck ist von einer bemerkenswerten Hinfachheit. Etwas unhandlicher ist der 
Potentialausdruck (9), da im Integranden eine relativ koniplizierte arc tg-Funktion auftritt. 
Der Ausdruck (9) ist aber zur Lésung der vorgesehenen Randwertaufgabe nicht unmittelbar 

forderlich. 7 
i Bemerken wollen wir noch, daB man mit c—> © aus (13) den bekannten Ausdruck fiir den 
Abwind eines gebundenen Wirbels variabler Intensitaét mit seiner Schleppe freier Wirbel im 
Inkompressiblen erhalt. 


e« 4. Der harmonisch pulsierende gebundene Wirbel mit seiner freien Wirbelschleppe im kom- 
pressiblen Medium. Wir machen fiir die Belegungsfunktion einen harmonischen Ansatz 
g(t—t)Hiadel-4 (14) 
27* 
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und erhalten damit aus (13): 
t+T 


K(x,1='@ eter f eine 8 1/) (")' ae. 


Wir denken uns hier T sehr groB. Das heif®t, wir wollen annehmen, das Wirbelsystem sei schon 
vor sehr langer Zeit (T'— co) entstanden. Wir betrachten zunichst den Fall x= — |x|. Dann 
liegt auf dem Integrationsweg kein Pol, und man kann schreiben 

f dE f -e6 dé), x= — [xfs 


ay t+ T 


Lo 


K (x, ) =" 


eit 


v 


Der zweite Term stellt fir T’—> 00 weiter nichts dar als den Abwind eines stationaren gebun- 
denen Wirbels. Er interessiert uns im Rahmen unserer Betrachtungen nicht. Wir kénnen ihn 
durch geeignete Anfangsbedingungen annulliert denken und lassen ihn deshalb in Zukunft immer 
weg. Es ist dann 


; SES: caaee Ie ae 
Ried eres hale fe yj eee fir x«x=— |x|, T+ o. ~ (15a) 


Wenn x= + |x| ist, liegt auf dem Integrationsweg immer der Pol = x/v, da im Falle v<e 
immer a)< x ist. Dieser Pol liefert einen besonderen Beitrag zum Integral. Man findet mit 
derselben Betrachtung wie vorher, und wenn man wieder den vom stationaren Wirbel herriihren- 
den Term weglaBbt, 


L = ae TE 
adh re a =a (emasiinlgy 
— \ 
= : (15b) 
awe 5) toe x= tel, To 


Die am Orte des gebundenen Wirbels (x= 0, y= 0) wirkende Auftriebskraft kénnen wir 
nach (14) und (10) sofort angeben. Wir subtrahieren verabredungsgem48 den stationaren Term 
nnd finden 
t 
A(t)=2z0ev[iw f eHtde ye ALOd re 2p vee’. (15c) 
7 
Wir versuchen jetzt, die Ausdriicke (15a) und (15b) auf Integraldarstellungen bekannter Funk- 
tionen zuriickzufiithren. Wir setzen dazu 
She] Mx 
 e (l= M3) (==) 


g mit M=—<1 


und fiithren noch die Bezeichnung ein 


K (x, t) = K (x) e'?, (16) 
Damit erhalt man 
a A Ls 3 io Mx i —— 1@ |x| u i@x 
——— (ED ‘ = ] eRe 
K ele ee c (1—m*) yut—i c(1I— M2) ae % +P Ed @ v 
() ae ies |x tie ie 2x ov . (17) 
Mx 


al 
Beim Wurzelausdruck ist der Waurzelzweig zu wahlen, fiir welchen gilt: 


jw l1—>-+u fir u>o. 
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Das Integral in (17) 1aBt sich weiter umformen. An Stelle von 


a a ad ' 
G= fe ae 7 eu mit k= + |k| 


il} Us 


kann man namlich schreiben 


t 00 nee : 
“G=— [ Vw—Teit f élvea) df du. 


Das Integral itber € existiert fiir alle Werte von u, fiir die z. B. gilt 
1 ] 
oe rye Re [u] < +i : 


Man kann in der Tat bei der Integration in der u-Ebene den Weg so von —1 nach i © fihren, 
das diese Bedingung auf dem ganzen Wege erfiillt ist (da M<1 ist). Dann kann man, wie sich 
zeigen laBt, die Reihenfolge der Integration vertauschen. Das Integral iiber u ist eine bekannte 
Integraldarstellung der Hankelfunktionen. Man findet 


== co 


c 
cat . MH (k —if) 
ares k—il 

0 


ae: 
Mit k— if = ka erhalt man 
_ tk Aer ka) pn ake 
= 45 tu Hee), M do, 
1 


wenn man beachtet, das man den Weg statt nach 1 + i © auch nach -} © fiihren darf. Man 
kann also an Stelle von (17) auch schreiben 


fo 2] 


tox 
1@MxX 1Mx 
= x inwe ” le w |x| 0 Saal A t=) CO eae ae 
— 2 2 vyi(1 M?) — % — a Be 18 
ae) x| ae Ht ry |° 3 cme pag cy) 
1 


Dieser Ausdruck 1aBt sich mit bekannten Formeln der Besselfunktionen weiter umformen, was 
fiir numerische Rechnungen vorteilhaft sein kann. Nach einigen Umformungen erhalt man 


schlieBlich 


O eMas. 7 
Sry ht TO tie | o [lal Me | o |x : } 
Ds oreer esa Tre | He qa | rei al [ed — M4) |f \ 
ioax 1 iM xX tw % (19) 
tm % Ww ae HP oo |x| o 20—M) dy tile * In———— ieee 
2 v2 V1 — M? te ie MP) 108 v 1—y1—™M? 
0 


Damit kennen wir die auf der x-Achse vorhandene in y-Richtung weisende Geschwindigkeits- 
komponente eines Wirbelsystems, welches im Koordinatenursprung den Auftrieb 


A= 2 mov exp (iwt) 
erzeugt. 
5. Die vom Medium der Schallausbreitung auf die harmonisch schwingende Platte ausgeiibten 
Kriafte. Wir denken uns jetzt eine Belegung der durch Gl. (19) beschriebenen Wirbelsysteme. Die 


Belegungsfunktion (x) sei nur im Bereich —1 S$ x S +1 von Null verschieden. Die dadurch 
auf der x-Achse erzeugte senkrechte Geschwindigkeitskomponente sei 


w (x,t) = we (aye! 


Dann ist 
Sie , 
w(x) = f K (x — %1) y (%) dx,. (20) 


—l 
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Der Belegungsstirke der Wirbel y(x) entspricht nach (15c) die Kraftbelegung 


p (xt) = p(x)! = 2apvy (a) ea (21) 
An Stelle von (20) kénnen wir also auch schreiben 
; +1 
10x) = ga | K (w— x) p(™) day. (22) 
=I 


Das ist, wenn w vorgegeben ist, eine Integralgleichung zur Bestimmung der Druckbelegung p(x). 
Eine beidseitig unendlich lange Platte, deren Vorder- bzw. Hinterkante bei «= — 1 baw. 
x= +1 liegt, fiihre kleine harmonische Schwingungsbewegungen in y-Richtung aus, die durch 


f (x,t) =f(x) eet mit [f|<1, |f'|<1 
beschrieben seien. Dadurch wird dem Medium an der Plattenoberflache der Abwind auf- 
gezwungen 


Wt (x) eb! = iw f (x) +o<f(x) eres. (23) 


In einer linearisierten Theorie darf man nun annehmen, da dieser Abwind an der Plattenober- 
flache mit dem Abwind auf der x-Achse (— 1 < « < +1) identisch ist. Damit ist w(x) durch 
die Plattenschwingung eindeutig festgelegt, und die auf der Platte wirkende Auftriebsverteilung 
kann durch die Integralgleichung (22) bestimmt werden. Die Eindeutigkeit der Lésung wird 
durch die Kuttasche AbfluBbedingung [p(-++1)= 0] gewahrleistet. 

Man iiberzeugt sich leicht, daB (22) mit der von Possio! fiir dasselbe Problem angegebenen 
Integralgleichung identisch ist. Damit haben wir mit Hilfe des neueingefiihrten Wirbelpotentials 
der kompressiblen Unterschallstrémung dasselbe Ergebnis erzielt, das Possio ausgehend vom 
Prandtlschen Beschleunigungspotential gefunden hat. Von der Integralgleichung (22) kennt 
man auer fiir den Sonderfall c— oo noch nicht die strenge Lésung. 


6. Zusammenfassung. 


Es wird gezeigt, daB die bei Strémungsproblemen des inkompressiblen Mediums so erfolgreich 
angewandte Wirbelmethode auf das kompressible Medium zur Untersuchung von ebenen Schall- 
ausbreitungs- und Strémungsvorgangen iibertragbar ist. Dazu werden fiir die Unterschall- 
strémung Geschwindigkeitspotentiale gewisser elementarer Wirbelfelder hergeleitet. 

Als erstes wird das Potential eines Wirbelsystems angegeben, das aus einem gebundenen 
Wirbelzentrum konstanter Auftriebskraft besteht und dem zugehérigen freien Wirbelzentrum, 
welches mit der Strémungsgeschwindigkeit des Mediums abschwimmt. Im Gegensatz zum 
Sonderfall des inkompressiblen Mediums 1a8t sich jedoch der Potentialausdruck nicht in zwei 
Summanden aufspalten, von denen der eine den gebundenen Wirbel allein und der andere den 
freien Wirbel (Anfahrwirbel) allein darstellt. 

Als nachstes wird das Geschwindigkeitspotential eines gebundenen Wirbelzentrums variabler 
Auftriebskraft angegeben mit seiner auf der Spur des gebundenen Wirbels liegenden Belegung 
freier Wirbel. 

Als drittes wird ein Sonderfall des zweiten Wirbelsystems untersucht, bei dem die im ge- 
bundenen Wirbelzentrum angreifende Auftriebskraft harmonisch mit der Zeit variiert. 

Zuletzt wird die praktische Brauchbarkeit dieser ,,kompressiblen‘‘ Wirbelpotentiale erwiesen, 
indem mit ihrer Hilfe die Integralgleichung fiir die Druckverteilung einer angestrémten, schwin- 
genden Platte unendlicher Spannweite und endlicher Tiefe hergeleitet wird. 


(Eingegangen am 20. April 1950.) 


Anschrift des Verfassers: Dr.-Ing. Johannes Dérr, Darmstadt, Steinbergweg 39. 


1 C. Possio, Aerotecnica 18 (1938), S. 441. 
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Beitrige zur Theorie der durch gleichférmigen Schub 
beanspruchten Platte. 


I. Mitteilung. 
Von W. Kucharski. 


1, Allgemeine Einleitung; kurze Inhaltsangabe. In der folgenden Mitteilung und in den 
weiteren Mitteilungen II und III zum gleichen Thema wird die bekannte Methode der so- 
genannten ,,Stérungsrechnung“ auf die lineare partielle Differentialgleichung einer durch gleich- 
formigen Schub beanspruchten Platte angewandt, bei welcher die Auslenkungen quer zur 
Mittelflache klein gegeniiber der Plattendicke sind, und diese wiederum klein gegeniiber den 
Langen- und Breitenabmessungen der Platte. Abschnitt 2 enthalt die notwendigen Einzelheiten 
tiber die allgemeine Problemstellung, die zu lésende Differentialgleichung, die benutzten Be- 
zeichnungen u.dgl. In Abschnitt 3 wird die gleichférmig verteilte Schubspannung auch als 
zeitlich konstant vorausgesetzt und das so entstehende Problem durch geeignete Ansitze geldst. 


Die allgemeinere Aufgabe, welche bei einer periodisch um einen konstanten Wert schwan- 
kenden Schubbeanspruchung entsteht, wird in Mitteilung II behandelt werden; dabei kommen 
nach geeigneter Vorbereitung fiir die hier vorliegende partielle Differentialgleichung vierter 
Ordnung einige Gedankenginge und SchluBweisen sinngema® zur Anwendung, die zur Lésung 
der analogen gewéhnlichen Differentialgleichung fiir ein System von einem Freiheitsgrad mit 
-periodisch veranderlicher elastischer Kraft gefiihrt haben. Hieran werden sich einige Ausfih- 
rungen vergleichender Art beziiglich der verschiedenen Liésungsansatze schlieBen. 


In Mitteilung III soll schlieBlich die grundlegende Differentialgleichung als Wellengleichung 
behandelt werden, woraus sich einige verhaltnismaBig einfache Lésungen fiir bisher nicht be- 
handelte Plattenformen und Gesichtspunkte fiir allgemeinere Aufgaben ergeben. 


Mit Ausnahme dieses letzten Abschnittes beziehen sich die Untersuchungen auf rechteckige 
Platten, an deren AuSenkanten die der jeweiligen Aufgabenstellung entsprechenden Schub- 
spannungen als dufere eingepragte Krafte angreifen. Die hier entwickelten Methoden lassen 
sich ohne wesentliche Anderungen auch auf diejenigen Falle iibertragen, in welchen an den 
Plattenseiten auBer den Schubspannungen (oder ohne diese) gleichférmig verteilte Zug- und 
Druckspannungen angreifen. 


Zur naheren Befassung mit den Problemen der Plattentheorie wurde der Verfasser u. a: 
durch die Bemiihungen von Herrn Dipl.-Ing. Thomas angeregt, der fiir das in Abschnitt 3 be- 
handelte Eigenwertproblem einige Rechnungen nach dem Ritzschen Verfahren anstellte. Herrn 
Dr. Metzmeier dankt der Verfasser fiir weitere Anregungen und wertvolle Literaturhinweise, 
besonders auf die interessante Arbeit von Herrn Dr.-Ing. Moheit, ,,Schubbelastung rechteckiger 
Platten“ (Leipzig 1939), welcher durch Differenzenrechnung und Versuche im groBen Mafstab 
iibereinstimmende Eigenwerte und Plattenformen feststellte, wahrend in frither verdffentlichten 
Arbeiten zum gleichen Problem aufgrund des Ritzschen Verfahrens zwar der Kigenwert ohne 
wesentlichen Unierschied, aber die zugehérige Plattenform mit grundsatzlicher Abweichung 
(eine anstatt zweier Beulen) angegeben war. Die zitierte Arbeit enthalt auch ausfiihrliche Lite- 
raturhinweise. 

Beziiglich der mathematischen Grundlagen sei z. B. auf die bekannten Biicher von Courant- 
Hilbert und Frank-Mises verwiesen; bei der Behandlung des auf eine partielle Differential- 
gleichung ausgedehnten Mathieuwschen Problems hat der Verfasser die Kapitel iiber die ent- 
sprechenden Fragestellungen in ,,Horn, Gewohnliche Differentialgleichungen“‘ benutzt. 

Das Ziel der vorliegenden Arbeiten des Verfassers besteht in der Herausarbeitung allgemeiner 
Methoden und Gedankengange fiir die unmittelbare Lésung der in der Literatur bisher wenig 
gepflegten partiellen Differentialgleichung der Plattenbewegung ; Einzelfalle werden hier ledig- 
lich als Beispiele fiir die Methode behandelt. Weitergehende numerische Berechnungen fiir die 
unmittelbare Anwendung befinden sich im Institut des Verfassers in Arbeit; hieriitber soll in 
besonderen Veréffentlichungen zu gegebener Zeit berichtet werden. 
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2. Problemstellung; Bezeichnungen; die partielle Differentialgleichung. Die rechteckige Platte 
von der Gesamtdicke h mit den Seitenlangen a und 6 sei nach Abb. 1 in dem rechtwinkligen 
(x, y)-System orientiert; die Auslenkungen w ihrer Mittelebene werden in der z-Richtung als 
positiv gerechnet, die mit den angenommenen x- und y-Richtungen ein rechtshandiges System 
ergibt. Es sei @ die spezifische Masse pro Volumeneinheit, t die Zeit und t die an den Randern 
in den gezeichneten Richtungen in gleichférmiger Verteilung angreifenden Schubspannungen (auf 
die Flaicheneinheit bezogen). Damit ist, wenn 7 nicht mit t variiert, ein gleichformiger Spannungs- 
zustand Tyy==Ty,—=7 in der gesamten Platte verbunden; das gleiche trifft mit geniigender 
Genauigkeit zu, wenn 7 sich mit t geniigend langsam andert, was im folgenden stets vorausgesetzt 
wird. Ferner sei h geniigend klein gegeniiber a und b, w geniigend klein gegeniiber h, so daB 
die bekannte lineare Plattentheorie gilt. Bezeichnet dann noch v= 1/m die Ziffer der Quer- 
kontraktion, E den Elastizititsmodul, und werden partielle Differentiationen nach den Ver- 
anderlichen x, y, tin der itblichen Weise durch Indices angedeutet, so lautet die lineare partielle 
Differentialgleichung fiir die Auslenkung w bekanntlich 

Eh? 
O Wey Ti (19) AAw 2 T Wz = 0 (1) 
mit 
NAW Wig te Wingy i Wenn 


Fiihrt man neue, dimensionslose Variable ein durch 


(2) 


| 


Yay Sarl, shee. wi as | 
| 


SANS AF ra RES 
eee, (ers) dato anal 
und schreibt dann zur Vereinfachung wieder w anstatt w,, x anstatt &, y 
Z anstatt 7 und tanstatt t,,so erhalt man mit 
__ 24(1—r?) @? 
am WE” (3) 


Abb. 1. oD 


anstatt (1) die grundlegende Differentialgleichung in der Form 
wy, +AAw-+0w,,= 0. (4) 


Hierin kann # eine beliebige Funktion der Zeit sein, mit der einzigen Einschrankung, daB 
dd/dt geniigend klein ist, daB man also die Massenbeschleunigungen tangential zur Mittelflache 
der Platte mit geniigender Genauigkeit vernachlassigen kann, wie das im folgenden durchweg 
geschehen soll. 

Gleichung (4) ist unter Beriicksichtigung der jeweils gegebenen Rand- und Anfangshedin- 
gungen bei gegebener Zeitfunktion # zu lésen. Als Randbedingungen kommen fiir die Anwendung 
vor allem in Frage: 


Ebene Auflage, d. h. w = 0 auf dem Rande fiir alle t (5) 
und entweder 


vollkommene Einspannung, d.h. w, == w,—= 0 auf dem Rande fiir alle t (6) 
oder 


momentenfreie Auflagerung, d. h. dw — 0 auf dem Rande fiir alle t. (7) 


Als Anfangsbedingungen, z. B. fiir t= 0, kommen, wie fiir jedes kinetische Problem, zwei 
unabhangige Bedingungen in Frage, z. B. 


w(x, ¥, 0) = W(x, y), | 


W(X V5 0) —= Ve (aaaty | 


d. h. Lage und Normalgeschwindigkeiten der Plattenpunkte fiir t— 0, gegeben durch die vor- 
geschriebenen Koordinatenfunktionen W bzw. V. 


und 


(8) 
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3. Die Schubspannung T ist zeitlich konstant. a) Aufspaltung der partiellen Differen- 


Melsioieh uns: In diesem Abschnitt wird 7 als zeitlich unveranderlich angenommen; in (4) ist 
jetzt also 


He konst, <0. (9) 
Es werde sofort die Frage nach stehenden Schwingungen gestellt durch den Ansatz 
tp = T, (1) Ug(% y) (0 = 0,1, 2, 3,-°+), (10) 
der in bekannter Weise zu den Differentialgleichungen fiir U und T fiihrt, niamlich 
AA U, +0 U pay — Ap Up = 0 (11) 
und 
Tai Ag do. 0 > (12) 


Gleichung (11) stellt bei gegebenem # und vorgeschriebenen Randbedingungen ein Eigen- 
wertproblem dar mit den Eigenwerten 4, und den dazugehérigen Eigenfunktionen U,. Ist dieses 
gelést, so liefert (12) sofort die Lésung 


T, = A, sin w,t + B, cos wot | 


mit (13) 
Op hos | 
Eine vollstandige Lésung hat hiernach die Form 
w= 5) (Ag sin wot + By cos wot) Up, (14) 


wobei die Integrationskonstanten A, und B, aus (8) zu bestimmen sind, nachdem W und V 
nach Eigenfunktionen U, des betreffenden Falles entwickelt sind. 

b) Zur Existenz der Lisung. DaB Gleichung (11) fiir die physikalisch in Frage kommen- 
den Randbedingungen lésbar ist, laBt sich in einer hier geniigenden Weise erharten wie folgt: 
Zweifellos hat 

AAU + 9U xy = p(*, 9) (15) 
eine Léisung, worin p einen beliebig verteilten Normaldruck auf der ruhenden Platte bedeutet. 


Wird hierbei p so gewahlt, daB es auBerhalb einer passenden Umgebungsflache F eines Punktes 


F 

Xq: ¥o Null ist, und wird dann diese Flache mit p+ 0 unbegrenzt verkleinert, wahrend [pat 
einen beschrankten, von Null verschiedenen Wert behalt, so ergibt sich in bekannter Weise die 
Belastung durch eine Einzelkraft im Punkte x,, yy. Hierfiir existiert also ebenfalls die Lésung. 
Diese ist linear; die Lésungen fiir beliebig viele Kinzelkrafte entstehen somit durch Superposition 
und sind unabhangig von der Reihenfolge, in welcher die Einzelkrafte aufgebracht werden. Das 
gleiche gilt fiir die potentielle Energie der Platte. Es gelten also die nach Maxwell usw. benannten 
Reziprozitatssatze. Das hei®Bt aber: Zu dem Differentialausdruck 


(y= AAY 8 Uy 


existiert eine symmetrische Greensche Funktion. Nach bekannten Satzen aus der 
Theorie der Integralgleichungen ist damit die Lésbarkeit von (11) sichergestellt. Man kénnte 
daher auch (11) mit Hilfe der zugehérigen Integralgleichung lésen, wenn es gelange, die Green- 
sche Funktion in einer brauchbaren Form aufzustellen. Einige Vorversuche haben jedoch gezeigt, 
daB dies auch im vorliegenden Falle nicht gerade zu einer Vereinfachung des Problems fihrt. 

c) Sonderfalle, die hier nicht in Frage kommen. Herr Dipl.-Ing. Thomas hat 
bemerkt, daB Gleichung (11)in ge wissen, vonihm nicht weiter untersuchten Fallen Lésungen 
von besonders einfacher Art haben kénnte. Schreibt man die Gleichung, in welcher iden- 
tisch verschwindend U ohne Interesse ist, in der Form 


AAs Oe a 
a eete a= 9" (11a) 
so werden Méglichkeiten nahegelegt, bei denen die Quotienten 
AAU air: vi 
a= <a und b = U ( ) 
yon den Koordinaten und von @ unabhangig sind. In solchen Fallen ware dann der Eigenwert 
A=a+bb, 


also linear in ?. 
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Es ist leicht zu zeigen, daB diese Méglichkeiten nur unter ganz besonderen, stark einschranken- 
den Voraussetzungen bestehen kénnen, die gerade in den wichtigsten Fallen der Anwendung 
nicht erfiillt sind. Durch (16) und die Voraussetzung a= 0 und by= 0 zerfallt namlich (11) 
in zwei gleichzeitig zu erfiillende Differentialgleichungen fiir U, von denen z. B. die eine 


Uy bU=0 j (17) 


lautet. Diese Gleichung mu durch U unter Erfiillung der Randbedingungen befriedigt werden. 
Wichtig ist vor allem die Randbedingung (5), wonach also die Platte aufihrem Rande festgehalten 
wird. Gleichung (17) ist hyperbolisch; sie tritt auch in anderen Problemen der angewandten 
Mechanik auf und ist in jeder praktisch wichtigen Richtung vollstandig erforscht, so daB die 
Lisung fiir die Randbedingungen (5) sofort hingeschrieben werden kann. Die Plattenrander, 
auf denen nach (5) der Wert von U mit Null vorgeschrieben ist, sind nun gerade die Charakteri- 
stiken der Differentialgleichung (17), und die bekannte Theorie dieser Gleichung liefert das Er- 
gebnis, da® U mit der Randbedingung (5) im gesamten Gebiet der Platte identisch verschwindet. 

In den praktisch wichtigen Fallen der Randbedingung (5) enthalten also die Quotienten a 
und b nach (16) noch ? und im allgemeinen auch die Koordinaten, und von einer linearen Ab- 
hangigkeit des Eigenwertes A von # kann hierbei nicht die Rede sein. 

Ob es physikalisch interessante Randbedingungen gibt, bei denen (17) zu erfiillen ist, bleibe 
hier dahingestellt. 

Zur Vermeidung von MiBverstandnissen sei noch folgendes vermerkt: Im allgemeinen wird /, 
als Funktion von # durch eine gekriimmte Kurve dargestellt werden. Diese 1aBt sich fiir eine 
geniigend kleine Umgebung eines Punktes ,,immer“‘ mit beliebiger Genauigkeit durch ihre Tan- 
gente ersetzen, etwa mit der Gleichung 1 = a, + 6,0. Die hierbei auftretenden ,,Konstanten™ a, 
und 6, haben aber nichts mit denjenigen nach Gleichung (16) zu tun und sind von Punkt zu 
Punkt der Kurve A = A(@) verschieden, also von # abhangig. Der Ersatz dieser Kurve durch 
ihre Tangente fiir die Umgebung des Berithrungspunktes bedeutet also nicht, da U in dieser 
Umgebung angenahert die Differentialgleichung (15) erfiillt. 

d) Der ,,Stérungsansatz“;eine besondere Schwierigkeit;zur Konvergenzfrage. 
Zu brauchbaren Lésungen von (11) gelangt man in an sich bekannter Weise durch einen so- 
genannten ,,Stérungsansatz beziiglich des Parameters ? von der Form 


he = log + ]oP + p08? + Ip98% + °°, | 

Ug = Ugo + Ug, P + Upe D” + Ugg? ++ -- 

In diesen Reihenansatzen nach Potenzen von @ sind die 1 Konstante, die u, Fuktionen von 
xund y. Geht man mit (18) in die zu lésende Gleichung (11), ordnet nach Potenzen von # und 


setzt den Faktor jeder Potenz von # fiir sich Null, so entstehen zunachst rein formal die folgen- 
den partiellen Differentialgleichungen, bei welchen der Index 9 zur Vereinfachung fortgelassen ist : 


(18) 


AA u,—Il, uy, = 0, 

AA wu, —1, uy — 1, uy = — Ugey 

AA u,— 1, u, — 1, ty = — Uyey +1, uy, 

AA u,—I1,u,—1, uy = — Ugzy +1,u, +1,u,, | em) 
AA wu —I1,u- —1, u=— uz.) +1,u, -+1,u, +1, 4, 


Diese unendlich vielen Gleichungen sind sukzessive zu lésen, wobei jede der Funktionen u 
fiir sich die jeweiligen Randbedingungen zu erfiillen hat. 

Die erste Gleichung fiihrt auf ein bekanntes Problem: J, bzw. u, ist offenbar ein Eigenwert 
bzw. die dazugehérige Eigenfunktion fiir die schubfreie Platte mit — 0. Diese Lésungen bei 
den Randbedingungen (5) mit (6) oder (7) sind bekannt und daher ohne weiteres zu benutzen. 

Die weiteren Gleichungen sind inhomogen, mit jeweils bekannten rechten Seiten, wahrend 
z. B. auf der linken Seite der zweiten Gleichung die Koordinatenfunktion wu, und die Konstante i 
unbekannt und zu ermitteln sind. Dies gelingt in der Weise, da jeweils die rechte Seite ebenso 
wie die unbekannte Funktion auf der linken nach den Kigenfunktionen der schubfreien Platte 
entwickelt wird, wobei die zunachst unbekannten Koeffizienten auf der linken Seite sich aus den be- 
kannten der rechten Seite Glied fitr Glied ergeben. Hierbei bleibt im allgemeinen auf der rechten 
Seite ein Glied mit der Eigenfunktion u, iibrig, das dann zur Bestimmung der noch unbekannten 
Konstanten I fiihrt. Diese ist Null, wenn ein solches Glied auf der rechten Seite nicht erscheint. 
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Im vorliegenden Falle tritt hierbei noch eine besondere Schwierigkeit auf: Es sei z. B. ein 
bestimmtes uy zugrunde gelegt, welches die vorgeschriebenen Randbedingungen erfiillt; u, wird 
dann als Reihe in Kigenfunktionen angesetzt, die ebenfalls die Randbedingungen erfiillen, wo- 
mit dann die gesamte linke Seite der zweiten Gleichung (19) den Randbedingungen geniigt. Der 
rechts auftretende Ausdruck uy,, wird aber so, wie er nach Differentiation von U, entsteht, im 
allgemeinen die Randbedingungen nicht erfiillen. Der Ausweg aus dieser Schwierigkeit besteht 
darin, da} das zunachst durch unmittelbare Differentiation entstandene Upxrx ebenfalls erst nach 
Kigenfunktionen, welche fiir sich die Randbedingungen erfiillen, entwickelt und dann erst rechts 
eingesetzt wird. Dies kénnte zundchst paradox erscheinen und fiihrt auch gegebenenfalls zu 
Reihen, welche in einem den Rand enthaltenden Gebiet nicht gleichmabig konvergieren. Diese 
und andere Einzelheiten werden spiiter an Hand eines Beispieles genauer erértert. 

Vorher sei kurz die Frage nach der Konvergenz des Verfahrens und der damit entstehenden 
Reihen behandelt. Da die Lésung fiir = 0 existiert und man — schon aus physikalisch plau- 
siblen Griinden — Stetigkeit an dieser Stelle in bezug auf # erwarten darf, wird mit Sicherheit 
ein von Null verschiedener Konvergenzbereich vorliegen. Die Frage nach seiner Erstreckung 
wiirde aber zur Beantwortung besonderer Untersuchungen bediirfen, die — soweit dem Ver- 
fasser bekannt ist — fiir derartige Probleme noch nicht allgemein vorliegen und bisher auch von 
ihm nicht in Angriff genommen wurden. Vom streng mathematischen Standpunkt aus bleibt 
daher eine wesentliche Liicke offen; aber es darf, wenn nicht zur Entschuldigung so doch zur 
Erlauterung, darauf hingewiesen werden, daB die Verhaltnisse bei manchen anderen Problemen 
der heutigen mathematischen Physik ganz ahnlich liegen, ohne daB man deswegen die so 
gewonnenen Liésungen als wertlos ansieht. 

Im folgenden wird, wie haufig bei derartigen Aufgaben, angenommen, daB ein gewisser 
Konvergenzbereich existiert, und simtliche Aussagen gelten unter der Voraussetzung, daf er 
nicht iiberschritten wird, wobei es zum Teil der praktischen Auswertung iiberlassen bleibt, die 
Grenze mehr oder weniger bestimmt herauszufinden. Bisherige numerische Auswertungen, die 
noch nicht abgeschlossen sind, lassen erkennen, dafi bei der quadratischen Platte und bei hier- 
von nicht zu stark abweichenden Seitenverhaltnissen der Konvergenzbereich itberraschend gro} 
zu sein und sorgar den kritischen Wert von #, der zum ersten Ausbeulen fithrt, wesentlich zu itber- 
treffen scheint, so daB auch noch negative Werte von A erfaBt werden. Dies ist physikalisch 
durchaus plausibel, wenn auch solche Lésungen nach kurzem Zeitverlauf ihren physikalischen 
Sinn verlieren. Genaueres hieriiber sei spateren Arbeiten vorbehalten. 

e) Ein Beispiel. Zur Erlauterung der Methode sei kurz die quadratische Platte mit 
den sogenannten Navierschen Randbedingungen (5 und 7) behandelt. Ihre Eigenfunk- 
tionen bei # = 0 sind bekanntlich 


Tua See Bit Vy 1A, Wa Loe OS es) (20) 
oder passende Kombinationen hiervon, mit den Eigenwerten 
Leen =e (21) 
bzw. lOpm,n(—— mens, 
Es sei fiir die Lésung von (19) von dem tiefsten Higenwert 
h=(l+1P=4 | Bs 
mit der Eigenfunktion Uy = sin “sin y | 
ausgegangen. Die erste der Gleichung (19) wird hierdurch erfiillt. Die zweite Gleichung lautet 
hiermit 
AA u,— 4u, — 1, sin x sin y = — cos cosy. (23) 
Durch die rechts stehende Koordinatenfunktion werden offenbar die Randbedingungen nicht 
erfiillt ; sie ist daher als Reihe in den Eigenfunktionen darzustellen. Diese Reihen sind bekannt ; 
es sei z. B. auf Knopp, Theorie und Anwendung der unendlichen Reihen, 3. Auflage, S. 390 ver- 


wiesen. Sie lauten 


cos px == 5 Fee. sinox(o = 1,3,5,...), p gerade, | 
ne | (24) 


oo 


e i 2.4.05" =) n de 
px=-— —— sin Ux (uU = 4, 4, O, , p ungera 
Pee a eer px 


IU 


“M=2 
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und gelten mit gleichmafiger Konvergenz im nicht abgeschlossenen Intervall 0 <x <a. Die 
Endpunkte sind auszunehmen, da die Funktion cos px durch die Entwicklung (24) als ungerade 
Funktion fortgesetzt wird und daher in den Endpunkten unstetig wird. Hierdurch entsteht 
fiir die Lésung keine Schwierigkeit. Denn die einzelnen Glieder der Reihenentwicklung fiir U 
werden so bestimmt, da® die Differentialgleichungen (19) auf dem Rande gliedweise erfiillt 
werden, und die Randbedingungen werden durch U auf dem Rande exakt erfiillt, da die ent- 
stehenden Reihen (im Konvergenzbereich beziiglich #) auch noch in den dritten Ableitungen 
gleichmaBig konvergieren. Hine gewisse Vorsicht ist nur bei Sonderfragen geboten, bei welchen 
z.B. AA U auf dem Rande erfordert wird. Die hierfiir durch viermalige Differentiation ent- 
stehenden Reihen konvergieren im abgeschlossenen Gebiet nicht mehr gleichmafig. Fiir solche 
wohl selten auftretenden Fragen miiSte man sich daher von innen dem Rande beliebig nahe an- 
nahern, um den stetigen Anschlu® an die Verhaltnisse im Innern nicht zu verlieren. 

Fiir die rechte Seite von (23) ist in beiden cos-Funktionen p= 1; es ist daher die zweite der 
Entwicklungen (24), d. h. diejenige nach geraden Vielfachen von x und y zu benutzen. Hiermit 
entsteht aus (23) 


AA u,— 4u,— 1, sin x sin y = BO, EST era een RY Lig V2 Ae Oo otc W opal oon 
be 


wobei also w und » die vollstandige Reihe der positiven ganzen Zahlen unabhangig voneinander 
zu durchlaufen haben. 
Fiir u, ist zunachst der allgemeine Ansatz 


Uy = Da Gna SI MBs NY (26) 


m, T™ 
zu machen. Hierin ist sicherlich 
C0) RE Bt ut — on, Pa elie 


da rechts in (25) keine Glieder mit ungeraden Vielfachen auftreten und die linke Seite, gleich Null 
gesetzt, héchstens mit n= m= 1] erfiillt werden kann. Die Gleichung 


AAu—4yu=—0 (27) 
hat namlich die einzige von Null verschiedene Lisung 
u= Asin xsin y (A beliebig) , (28) 


da 4 der hierzu gehérige Eigenwert ist. 
Setzt man daher 


u,=— Asin xsiny + 3), a,,sinuxsin vy , (29) 
Ms v 
so folgt aus (25) 
A beliebig , (30) 
1=0, (31) 
16 by 


Ft 8 [ue 9 A] (1) ae Co 


Das Ergebnis (30) bedeutet nichts anderes, als daB zu den entstehenden Liésungen der inhomo- 
genen Gleichung (29) noch beliebige Lésungen der homogenen Gleichung (27) hinzugefiigt werden 
kénnen. 
Aus (31) folgt, daB in der Potenzreihe fiir 4 nach (18) das lineare Glied mit % ausfallt; es 
ist also 
Of 
ao 
Die nach (32) entstehende Reihe fiir 4, konvergiert so gut, da® sie fiir Zahlenrechnungen 
ohne weiteres brauchbar ist; bei konstantem y gehen die Fourierkoeffizienten wie 1/u? nach Null 
bei u—> co. 
Fiir die dritte der Gleichungen (19) erhalt man aus (29) 


OF Mire? 07 (33) 


Ujxy = Acosxcos y + DS) uy ya,,c0s wx cosvy. 
uy 


Dies ist wieder nach den Eigenfunktionen unter Benutzung von (24) zu entwickeln und rechts 
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einzusetzen. Mit (31) geht dann die dritte Gleichung (19) iiber in 


A a2 ee : ; oe gl pw : 5 
A uy — 4, — 1, sin x sin y = 42 GOS sin 4% sin vy | 


16 o 34 
— >» (» uy ,a, pe @ (34) 
Bs. ( f 1 uy) = (02 — 2) (2 — 


(0,,G==1, 3. be: :) (CAD Re er 


psinexsinay | 


wobei in der letzten Doppelsumme die Reihenfolge der Summationen bereits vertauscht wurde. 
In dieser Doppelsumme ist mit 9 = o = | als erstes Glied dasjenige mit sin x sin y enthalten; 
links hat lediglich das letzte den gleichen Faktor; so erhalt man zuniichst 


16 yp 
L=2> ae 
2 8s (uw? =) G1) | 


(35) 


256 2 2 
c , Ler 
at 2 [(u® + 98)? — 4] (a? — 1) @? — 1) 


Von der Summe itber 0, cist hiermit das erste Glied verbraucht; es ergibt mit 1, sin x sin y auf 
der linken Seite Null. Das erste und zweite Glied links bildet mit der ersten Summe rechts 
bis auf den Faktor A die fiir u, geléste Gleichung (25) mit 1, = 0; u, enthalt also einen additiven 
Anteil Au,. Ferner ist die zweite Summe rechts zu beriicksichtigen, und schlieBlich ein beliebiger 
Anteil B sin x sin y fiir den homogenen Ausdruck der linken Seite zu addieren. So entsteht 


Uu,= B,) + Au, +» ,a,,sinoxsinoy | 
: O90 
mit a (36) 
wae Qo i pw 
ae (g" + 07)? —4 a? De 19%» (i — 92) (p® — 0?) (0,6 = 3,5, 7,°*"), | 
wy 


wobei der Strich am Summenzeichen an den Fortfall des Gliedes mit 9 = o = | erinnern soil. 
In sinngemaB gleicher Weise ist die Rechnung fiir die weiteren Gleichungen (19) durchzu- 
fiihren. Man iiberzeugt sich leicht, daB in diesem Falle stets 


0 (m ungerade); 


in den Potenzreihen nach # kommen also in diesem Spezialfall nur die geraden Potenzen vor. 

Wie man bereits aus (36) deutlich ersieht und fiir die weiteren Gleichungen leicht bestatigt, 
fiihrt die jeweilige Hinzufiigung von homogenen Lésungen mit unbestimmten Koeffizienten zu 
nichts anderem, als daB die Lésungen von (11) nur bis auf einen beliebigen Faktor bestimmt sind. 
Die Allgemeinheit wird also nicht eingeschrankt, wenn von vornherein 


A=B=0 (36) 


gesetzt und die allgemeine Vorschrift eingefiihrt wird, dai Lésungen der jewei- 
ligen homogenen Differentialausdriicke auf den linken Seiten der Gleichungen 
(19) zu den in der angegebenen Weise resultierenden Lésungen der vollstan- 
digen Gleichungen in keinem Falle hinzugefiigt werden sollen. 

Die Funktionen Uo nach Gleichung (11) sind damit fiir den hier behandelten Sonderfall voll- 
kommen bestimmt. 

Fir dasjenige Uo, das sich in der skizzierten Weise aus der Grundfunktion uy = sin » sin y 
ergibt, sind in Tabelle 1 die Werte der ersten Reihenkoeffizienten 4,» bzw. n@), angegeben, 
wie sie von Herrn Dipl.-Ing. Erlenwein in Vervollstandigung der ersten Ziffernrechnungen des 
Verfassers mit durchaus ertraglicher Rechenarbeit ermittelt wurden. Fiir den zugehérigen 
Eigenwert / ergibt sich daraus die Reihe 


= 4 — 9,146 - 10-3 92 — 5,28 - 10-* 04 — 3,57 - 10-28 s, (37) 
Fiir 1 = 0 ergibt sich hieraus der erste Wert von 0, bei welchem die Platte ausbeult, zu 
Drritisch — 18,71 . (38) 


Der bisher veréffentlichte, nach dem Ritzschen Verfahren ermittelte Wert ist eer 
Damit ist das Problem fiir die tiefste Eigenfrequenz der durch # verzerrten Platte voll- 


standig gelést. 
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Tabelle 1. yar (—10-). 


2 4 6 8 10 12 14 16 18 


2 1200,8437 | 72,7784 | 11,6085 2.9705 1,0096 0.4141 0.1940 0.1003 0,0560 
4 72,7784 | 11,3021 2,7448 0,8583 03246 0,1417 0,0691 0,0367 
6 11,6085 2,7448 0.9197 0,3530 0,1518 0,0720 0,0371 
8 2.9705 0,8583 0.3530 0.1596 0,0773 0,0399 
10 1,0096 | 0,3246 0.1518 0,0773 0,0414 
12 0.4141 0,1417 0.0720 0,0399 
14 0,1940 | 0,0691 0,0371 
16 
18 


1 — —150,354 | —11,596 | —2,066 | —0,572 | —0,206 | —0,088 | —0,04 
3 —150,354|+ 73,409} 110,454 | +2,326 | +0,714 | +0,272 | +0,122 | +0,060 
SD — 11,596 |+ 10,454] + 2,725 | +0,828 | +0.299 | +0,126 | +0,070 
7 — 2,066|+ 2,326) + 0,828 | +0,320 | +0,135 | +0,064 
9 — 0,572|+ 0,71 | + 0,28 +0,13 +0,06 
11 — 0,200)-- 0,28 | -- 0,13 +0,06 
13 
15) 


— 0,08 |+ 0,14 | + 0,06 
— 0,04 |+ 0,06 


2 —5,73 O59 t +0,11 +0,02 

4 +0,97 —0,16 —0,03 
6 +0,11 —0,03 
8 +0,02 


Auf Ergebnisse fiir die héheren Eigenwerte und die dazugehérigen Plattenformen sowie auf 
gewisse Eigentiimlichkeiten, die sich bei der erweiterten Rechenarbeit ergeben, soll spater in 
besonderen Arbeiten eingegangen werden. 

Hier sei noch darauf hingewiesen, daB mit den EKigenwerten eines Problems auch seine 
Greensche Funktion K(é, 7, x, y) gegeben ist, das bedeutet fiir die physikalische Anwendung die 
statische Auslenkung, die von einer konzentrierten, an der Stelle €,7 aufgebrachten Normal- 
kraft an der Stelle x, y hervorgerufen wird. Nach bekannten mathematischen Satzen ist 


Uo 5 Uo 3 
ee a (39) 


Auch dieses Problem ist durch die hier entwickelte Methode im Konvergenzbereich gelést. 

SchlieBlich sei vermerkt, dafi Uo durch den Ansatz (18) als Potenzreihe in # geliefert wurde, 
mit den Koordinatenfunktionen u,,, als Koeffizienten. Im Konvergenzbereich ist Umordnung 
erlaubt, z. B. nach sin m x sin n y , so dai Up dann als Fourierreihe in den Kigenfunktionen der 
unbelasteten Platte erscheint. Hierbei sind die Koeffizienten Potenzreihen in #. 
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f) In dem kurz behandelten Beispiel des vorhergehenden Abschnittes ergab sich dA/d# = 0 
fir = 0. Dies gilt aber nicht allgemein. So findet man z. B. nach der angegebenen Methode, 
daB durch den Punkt # = 0, = 1, = 25 drei Kurven A(#) hindurchgehen: 


7 Sy 64 be ; f ‘ 
a) eine mit (’(0) = — 9,2 fir up=sin xsinly +sinlxsin y, 


p ey, 64 cS : : : ‘ 
b) eine mit 4’(0) = Ir gga fiir uy = sin xsinly — sinlasiny, 


c) eine mit 4’(0) = 0 fiir uy—sinxsinly undfiir u y=sinlxsiny. 


Alle diese Méglichkeiten haben physikalische Bedeutung. Zur Veranschaulichung sei erwahnt, 
daB z. B. ug nach a) als Knotenlinie die Diagonale vom Punkte a, 0 nach dem Punkte 0, b hat, 
dagegen uy nach b) diejenige von 0,0 nach a,b. Der Ubergang von a auf b bedeutet also eine 
Drehung der Knotenlinie um 7/2. Dem Ubergang von # auf — # entspricht aber eine Drehung 
des Schubspannungspaares um 7/2. Geschieht beides gleichzeitig, was einem Ubergang von a) 
auf b) und von # auf — # entspricht, so muB 1() eine gerade Funktion, 4’ (0) == 0 demnach 
unstetig sein. Geht man dagegen innerhalb von a) oder von b) von +? auf — # iiber, so bleibt A’ (0) 
stetig. 

Da hiernach /4’(0) > 0 durchaus vorkommen kann, so wird die Platte nicht immer durch 
? > 0 im Sinne eines verkleinerten Eigenwertes ,,weicher“* gemacht, wie es bei u, = sin x sin y 
und ahnlichen einfachen Grundfunktionen der Fall ist. 


(Eingegangen am 8. Mai 1950.) 
Anschrift des Verfassers: Professor W. Kucharski, Berlin-Charlottenburg, Hardenbergstr. 34. 
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Beitrige zur Theorie der durch gleichformigen Schub 
beanspruchten Platte. 


II. Mitteilung. 
Von W. Kucharski. 


1. Periodische Schubbelastung; das Mathieusche Problem bei einer partiellen Differential- 
gleichung. a) Die zu lésende partielle Differentialgleichung. In der folgenden Mit- 
teilung II wird unter Hinweis auf die allgemeinen Problemstellungen und Grundlagen der im 
gleichen Heft dieser Zeitschrift veréffentlichten Mitteilung I zum gleichen Thema die Schub- 
spannung T als periodische Funktion der Zeit angenommen, die um einen konstanten Mittelwert 
pulsiert. Die technische Wichtigkeit dieser Problemstellung liegt auf der Hand. Mit 


0= 0 +0, sin wt (1) 
geht die grundlegende Differentialgleichung (4) von Mitteilung I iiber in 
wy, +A Aw +(O +) sin wt) w,y = 0. (2) 


Hierin sind bis auf weiteres 0,0, und w beliebige reelle Konstante mit der einzigen Kinschrankung 
beziiglich der Frequenz w, daB die in der Platte infolge der Schwankungen von Tt auftretenden 
tangentialen Massenbeschleunigungen mit geniigender Genauigkeit vernachlassigt werden kénnen. 
Eine Abschatzung hierfiir sei einer besonderen Untersuchung vorbehalten. 

b) Der ,,Stérungsansatz; Vorbereitungen zur Lisung der entstehenden Glei- 
chungen. Der ibliche Produkten-Ansatz nach Gleichung (10), Abschnitt 3 der Mitteilung I 
fiihrt hier nicht zum Ziele, wie man nach dem in Teil c) des genannten Abschnitts fiir eine ahn- 
liche Fragestellung durchgefiihrten Gedankengang leicht feststellt. Dagegen wird eine Lésung 
von (2) erméglicht durch den ,,Stérungsansatz‘ nach dem Parameter #, 


peas S 1083 = Wy + WD, + w,8§ + °°, (3) 


der nach Einsetzen in (2) und Zusammenfassung gleicher Potenzen von %, den folgenden Satz 
von unendlich vielen Differentialgleichungen liefert: 


wo, +Adm + Ow, = 0. (4) 
Wm, +A wn + O Wm, = = Wret o SUSAR (7 es Per) (5) 


Gleichung (4) mit O = konst. ist die in Abschnitt 3 der Mitteilung I geléste Differential- 
gleichung; man hat also sofort 


Cy Dy (A, sin Wot + B, cos wet) Uy (x, y), (6) 


o=0 
wobei sich U, aus der Differentialgleichung 


AAU, +00), —4yUp= 0 (7) 


ergibt und 2, =, der Eigenwert von (7) ist, zu welchem U, bei den vorgeschriebenen Rand- 
bedingungen als Eigenfunktion der Platte bei konstanter Schubbelastung O gehért. Hiernach 
sind im folgenden die U, und A, als bekannt anzusehen. 

Die Gleichungen (5) sind inhomogen; die rechte Seite wird jeweils durch die vorhergehende 
Gleichung als bekannte Funktion von x, y, t geliefert; sie sind sukzessive zu lésen. Wird Wo 
mit einem Faktor multipliziert, so multipliziert sich damit jedes w,, mit dem gleichen Faktor, 
und jedes w,, setzt sich linear und additiv aus allen denjenigen Einzelanteilen zusammen, die 
aus den U,sin@,t und den U, cos w,t fir 0 = 0,1, 2, ... flieBen, wobei jeder dieser Anteile 
fiir sich zu bestimmen ist, ohne Zusammenhang mit den zu einem anderen Index o gehérenden 


Anteilen. 
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Nun sei ,W,, diejenige Funktion, die aus dem wy-Anteil U, sin wt entsteht, 
ome ” ” ” rh oe) rr) uF cos Wot ” 


dann entsteht fiir w,, die Form 


Wm = SX (Ag* oWmy + Bo + omg)» (8) 
Q 


und fiir die vollstandige Lésung w erhalt man: 


w= So" DINE AR eS oe i Pec ms': (9) 
m a) m @ 


Hierin ist im Konvergenzbereich sicherlich die Umordnung erlaubt, die zu 


w= >) AS >) Wn OS, B, ome Dy 
@ m Q 


™m 


fiihrt, oder mit 


(10) 


zu der einfachen Form 
w= S (Awe + Bow, . (11) 
Q 

Hierin sind w,, und w,, Einheitsfunktionen, von denen w,, aus dem wy-Antell U, sin @,t 
und w,, aus dem w,-Anteil U, cos w,t entsteht. 

Dabei erfiillt jedes w,, bzw. wp, fiir sich die Differentialgleichung (2) und die von t unabhan- 
gingen Bedingungen auf dem Rande, die hier dieselben sein sollen wie in Abschnitt 3 der Mit- 
teilung I. 

Die Funktionen ,w,,, bzw. pwn. werden sich aus den inhomogenen Gleichungen (5) in ahn- 
licher Weise ergeben, wie die Funktionen u des vorigen Kapitels aus ihren inhomogenen Diffe- 
rentialgleichungen. Hierfiir wird festgesetzt, da®B Lisungen der betreffenden homo- 
genen Gleichungen in keinem Fall zu den sich bei dem Integrationsverfahren 
unmittelbar ergebenden Lésungen hinzugefiigt werden sollen. Dann gilt der Satz: 


Die Zeitfunktionen in ,w,,, und ,w,, enthalten keine anderen Frequenzen als @,und@. (12) 


c) Einzelheiten zur Herstellung der Lésungen. Die Liésung der Gleichungen (5) 
gestaltet sich im einzelnen wie folgt: Fiir m=1 erhalt man fiir ,w,, bzw. ,wj, die Gleichungen 


: é 1 
011, tA ow + OoWrr yy = — Uo,, sin Mgt sin wt—= —- > Uo,, cos (MW) — w) t | 
1 
+— U, cos (w,+w)t, | 
Wy, + AAW, + Oowys,, = — Up,, cos Wot sin wt = > Uo, Sim (Wp — w) t 


1 ; 
omer Uo,, sin (Wo @)t, 


zu deren Lésung offenbar mit Koordinatenfunktionen u, welche die Randbedingungen erfillen, 
anzusetzen ist: 
gy = gly 11 COS (pg — W) t + ollyy9 COS (pe +) Tt, \ 


oWy5 = ollyo1 SIN (Wy — W) t — Uyy2 SID (We --o)t, J 


(14) 


wobei sich die u nach Einsetzen von (14) in (13) in nunmehr bekannter Weise ergeben, und zwar 
eindeutig, da ein fiir allemal Lésungen der homogenen Gleichungen nicht hinzuzufiigen sind, so 
daB Satz (12) bestehen bleibt. Bei der Ausrechnung zeigt sich, daf je zwei der Funktionen in (14) 
einander gleich sind; es sind also nur zwei Koordinatenfunktionen zu bestimmen. 
Mit der so gewonnenen Lésung fiir m=—1 ist dann die rechte Seite fiir die Gleichung mit 
m=—2 zu bilden; diese selbst ist dann ahnlich wie soeben fiir m—1 skizziert zu behandeln. Bei 
der Umformung der trigonometrischen Zeitprodukte ergibt sich dabei die zunachst als stérend 
erscheinende Besonderheit, daB rechts Glieder mit sin w,t bzw. cos Wot auftreten, so da dann 
zusammen mit der linken Seite der sog. Resonanzfall vorzuliegen scheint. Der Ausweg aus 
28 
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dieser Schwierigkeit, die sonst das gesamte Verfahren sofort stoppen wirde, besteht einfach 
darin, daB die fiir diesen Fall bekannten stetigen Lésungen mit t cos wt bzw. t sin ot zu benutzen 
sind. Mit diesen geht dann das Verfahren gesetzmaBig weiter, wobei in den Gleichungen fiir 
m>2 naturgemaB Glieder mit #sinwt bzw. cos wt auftreten mit p>1 als ganze Zahlen. 

Das Auftreten dieser Potenzprodukte bedeutet nicht etwa, das die Lésung w mit wachsen- 
dem t auf jeden Fall unbeschrankt zunimmt. Es treten unendliche Reihen solcher Funktionen 
auf, iiber deren Verhalten zunachst nichts ausgesagt werden kann. Die Untersuchung dieser 
Frage wird naturgemaB eine der Hauptaufgaben dieser Mitteilung sein. An dieser Stelle des 
Gedankenganges ist festzuhalten, daB die Funktionen wo, und w,, in der skizzierten Weise ein- 
deutig und unter Erfiillung von (12) hergestellt werden kénnen und daher fiir das folgende als 
gegeben anzusehen sind. 

d) Vollstandigkeit der angegebenen Lésung. Ist man nun imstande, mit den ange- 
gebenen Funktionen die iiblichen Anfangsbedingungen eines kinetischen Problems zu lésen, 
mit anderen Worten: hat der vorgenommene Aufbau und die Durchfiihrung zu einer vollstan- 
digen Lésung gefiihrt ? 

Es seien z. B. Lage w(0) und Geschwindigkeit w, der Plattenpunkte fiir t—0 als Funktionen 
der Koordinaten gegeben. Diese und samtliche sonst auftretende Funktionen seien nach Eigen- 
funktionen U,(u= 0, 1,2...) entwickelt, mit den Koeffizienten 


fiir ane Moy), ? Meo, ? | 

w(0), We,(0), We,(0) 
und den Koeffizienten (15) 
i Trea Rt pea de | 

10,(0),  t,,(0)>  W,,(0) + | 


Alle diese Koeffizienten sind als bekannt anzusehen. 

Andererseits hat die Lésung w die Form (11) mit den A, und B, als Integrationskonstanten 
in der Anzahl 2 0. 

Setzt man die Entwicklungen mit den Koeffizienten nach (15) in (11) bzw. in den daraus 
gebildeten Ausdruck fiir w, ein, so erhalt man fiir die A, bzw. B, den Gleichungssatz 


SAieg, eS Been a | 
@ Q 
a Ao Boy, oe Bo Boy, = b, ? 
g 2 


Die Menge der y ist hier die gleiche wie diejenige der 0; (24) stellt also 2 9 lineare inhomogene 
Gleichungen fiir die 2 @ Unbekannten A, und B, dar. 

Nun ist zwar 9 unendlich. Man weil aber, da® derartige unendlich viele Gleichungen bei 
dem hier in Frage kommenden Aufbau der Koeffizienten einen von Null verschiedenen Konver- 
genzbereich besitzen, in welchem die Lisungen existieren und z. B. in bekannter Weise durch 
sukzessive Annaherung od. dg]. mit geniigender Genauigkeit bestimmt werden kénnen. 

Fir diesen Bereich ist dann die angegebene Lésung in dem hier in Frage (17) 
kommenden Sinn vollstandig. 

Es ist weiter physikalisch plausibel und wird ohne besonderen Beweis zu unterstellen sein, 
da8 die Liésung auch eindeutig ist. 

e) ,,Stabilitat der Bewegung. Bei der praktischen Anwendung wird man sich nur selten 
fiir die Einzelheiten der auftretenden Schwingungen interessieren. In den meisten Fallen wird 
die Hauptfrage diejenige sein, ob die durch eine Stérung eingeleitete Bewegung im Laufe der 
Zeit beschrankt bleibt (Stabilitat), wobei sie periodisch oder in sonstiger Weise mit konstant 
bleibenden Amplituden verlaufen oder auf Null abklingen kann, oder ob sie mit wachsendem t 
unbeschrankt anwachst. Diese Hauptfrage kann jetzt fiir die vorliegende partielle Differential- 
gleichung beantwortet werden. 

Dabei ist alles Bisherige derart formuliert, daB gewisse SchluBweisen, die fiir die Liésung der 
gewohnlichen Differentialgleichung 


2 +(a+bsinwt) z= 0 (18) 
bekanntermafen benutzt werden, auch sinngemaB auf die hier vorliegende partielle Gleichung 
angewandt werden kénnen. Das Wesentliche hierfiir war die von (9) auf (11) fiihrende Umord- 
pee verbunden mit der Einftthrung der Einheitsfunktionen Wp, und wy, welche den Satz (12) 
erfiillen. 


(ee eee ey: (16) 
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Nimmt man hierzu die Feststellung (17), so ersieht man, da® von den Einheitsfunktionen Wo, 
und wp, fiir ein bestimmtes @ zu einer Lisung auch alles beigetragen wird, was in den Zeitfunk- 
tionen die Frequenz @, und deren ganzzahlige Vielfache enthalt, genauer gesagt, da alle Anteile, 


die von den aus U, sin w,t baw. U, cos wot flieBenden Funktionen zu einer Liésung beigesteuert 
werden kénnen, in der Form 


aa We, + Bow, (19) 


enthalten sind, wobei A, und B, unabhiangig voneinander die reellen Zahlen durchlaufen, daB 
hierin keine anderen @, in den Zeitfunktionen auftreten, und da das jeweilige @, in den Zeit- 
funktionen mit anderem g nicht vorkommt. Zur Beantwortung der Stabilitatsfrage ist also jedes 
Paar der Einheitsfunktionen mit einem jeweils bestimmten 9 fiir sich zu untersuchen, und jedes 
Paar liefert zur Gesamtlésung einen Beitrag von der Form (19). Hierauf aufbauend ergibt sich 
folgender Gedankengang: 

Die grundlegende Differentialgleichung (12) bleibt ungeiindert, wenn man in ihr t durch 
t-+-P ersetzt, worin P= 2a/q@ die primitive Periode von sinwt bedeutet. Ist daher w eine 
Lésung von (2), so ist es auch w, wenn w aus w dadurch hervorgeht, daB hierin t durch t--P 
ersetzt wird. Infolgedessen sind auch diejenigen Funktionen w,, und w,, Lésungen, die aus We, 
bzw. wp, beim Ersatz von t durch t + P hervorgehen. 

Durch die Operation: Ersatz von t durch t-++P werden @, und seine Vielfache, die in Wo, 
und w,, vorkommen, nicht geandert, insbesondere nicht in entsprechenden GréSen mit anderem 
Index tibergefiihrt. Das heift aber: wy, und wo, gehdren ebenfalls zu der durch (19) dargestellten 
Mannigfaltigkeit von Lésungen; sie entstehen aus we, und wo, durch lineare Substitution, 


so daf 


Wo, = 411 Wo, + 442 Wo, » \ 


= (20) 
0) = 4q1 Wo, +499 Wo, J 
Hierin sind die Koeffizienten Konstante, deren Determinante nicht verschwindet: 
Gnd 
11 “12 
A= +0, (21) 
451 U0 


da die Gleichungen (20) naturgem4B nach w,, und we, auflésbar sind. Ferner sind die Koeffi- 
zienten als bekannt anzusehen; denn man hat lediglich die Lésungen w,, und w,, in der ange- 
gebenen Weise zu bestimmen, darin t durch t+ P zu ersetzen und die Substitutionskoeffizienten 
durch Vergleich festzustellen. 

Von solchen Substitutionen ist folgendes bekannt: Aus den urspriinglich gegebenen Grund- 
funktionen w,, und w,, lassen sich durch eine ganz bestimmte lineare Substitution andere Grund- 
funktionen Wo, bzw. Wo, bilden, derart, daB dann z. B. 


= sw, worin s= konst. (22) 


—_—; 
Ww 01? 


Q1 


Diese Konstante ist aus den Koeffizienten einer beliebigen Substitution aus deren sog. Funda- 
mentalgleichung zu bestimmen, welche lautet 


s? — (dy; +g) $ + (441 %2 — 49451) = 0. (23) 


Damit ist das Kriterium fiir die Stabilitat der Bewegung gefunden; denn 


die Bewegung klingt ab bei 0<|s|<1, wenn s reell; 

ex wachst an bei sia wenn s reell; 

a i ist periodisch bei s=-+1, mit der Periode P; (24) 
x - ist periodisch bei c= — 1s mit der Periode 2 P; 

i “A bleibt beschrankt bei |s|= 1, wenn s komplex. J 


Dieses Kriterium laBt sich einfacher formulieren. Weiter unten wird gezeigt, das bei der 
hier vorliegenden Differentialgleichung (2), in welcher das Glied mit w, fehlt, stets 


(25) 


041423 — 44249, = I. 
Hiermit lassen sich die Wurzeln von (23) auf die Form bringen 
— ex i2ah (26) 
28* 


Si 
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mit : 
cos 2a7h= im: (a,, +45) ; (27) 
und das Kriterium fiir Stabilitat lautet einfach 
1 
G} (44; +4)/ <1. (28) 


Durch (28) werden bei passender Diagrammdarstellung Kurven bestimmt, welche die Gebiete 
der Stabilitat von denen der Labilitat trennen, wie sie z. B. fiir die gewéhnliche Differential- 
gleichung (18) in der Physik und Technik als Struttsche Karte bekannt sind. Bei der hier vor- 
liegenden partiellen Differentialgleichung gibt es fiir jedes w, eine solche Karte, im ganzen also 
einen Atlas mit unendlich vielen Blattern. Weiter unten wird eine solche Karte fiir einen ein- 
fachen Fall entwickelt werden. 

f) Zwei Hilfssatze. Die praktische Berechnung der Einzelheiten fiir das Kriterium wird 
durch zwei Hilfssatze erleichtert, von denen der eine bereits erwahnt wurde. 

1. Die Ermittlung der Grundfunktionen we, und we, ergibt, daB sie sich in der Form 


Wey Ka (p= 0 es | 
P 


(29) 
wor = 5) K, Tp, | 
P 
schreiben lassen, worin die K, linear unabhangige Koordinatenfunktionen sind, die in beiden 
Reihen an der gleichen Stelle auftreten, wahrend die Zeitfunktionen is verschieden sind. 
Beim Ersatz von t durch t-+-P werden die Koordinatenfunktionen nicht geandert; dabei 
entstehen lediglich neue Zeitfunktionen, so daB 


Wo, i > K, Tp, | 
” f | (30) 
pe I ke | 

P 


Andererseits gilt die lineare Substitution (20). Da die K, linear unabhangig sind, kénnen 
aus (30), (20) und (29) keine Bedingungsgleichungen fiir die K, resultieren; also wird fiir jedes p 


Ty, = 441 Tp, +442 Tp, = Thy» \ (31) 
Th, = 41 T>, + Go T,, = dA es J 


das hei®t aber: 


Beim Ersatz von t durch t + P erfahren samtliche Zeitfunktionen in (29), und zwar jede | 
fiir sich, die gleiche lineare Substitution, wie sie die gesamten Funktionen w,, und wo, | (32) 
erfahren. | 

Fiir die praktische Ausrechnung der Substitutionskoeffizienten erbringt dieser Hifssatz den 
Vorteil, das man lediglich diejenigen zwei Zeitfunktionen zu verfolgen hat, die als gemeinsamen 
Faktor eine passend gewahlte Koordinatenfunktion besitzen. Da das Lésungsverfahren mit U, 
beginnt, wird man dieses als K, benutzen. Bei einem bestimmten Rechenaufwand erhalt man 
auf diese Weise beim Abbrechen der Rechnung die gréfte bis dahin mégliche Vollstandigkeit 
fiir die Zeitfunktionen und damit fiir die Substitutions-Koeffizienten. 

2. DaB die Determinante in Substitutionskoeffizienten allgemein gleich Eins 
ist, ergibt sich folgendermaBen: Die Grundlésungen w, und we, stellen fiir jedes konstant ge- 
haltene Wertepaar x, y, d.h. fiir jeden Punkt der Plattenmittelflache, zwei Funktionen der 
Zeit dar, die stetige Ableitungen erster und zweiter Ordnung nach der Zeit besitzen und von- 
einander linear unabhangig sind. Man kann von diesen Funktionen fir jedes konstant gehaltene 
Wertepaar x, y die Wronskische Determinante bilden, die wegen der Unabhangigkeit der Funk- 
tionen nicht verschwindet: 

We, We, 


W (wWo,, We.) = =~ 0. (33) 


Wor, w,, 


DaB wo, und w,, auBer von t auch noch von x, y abhangen, spielt fiir das folgende keine Rolle; 


man kann z. B. x, y als Parameter betrachten, denen man verschiedene Werte beilegen kann, 
die aber konstant zu halten sind, wenn es allein auf die Abhangigkeit von der Zeit ankommt. 
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Nach einem bekannten Satz aus der Theorie gewohnlicher Differentialgleichungen (1) gibt 
es nun eine und nur eine homogene lineare Differentialgleichung zweiter Ord- 
nung, fiir welche zwei derartige Funktionen ein Hauptsystem bilden, namlicn 


W (w, Wo1> Woo) 0 = ) 


W (Wo, Wo.) Oe 
mit 
w Wo, We, | 
W (w, Wo W,) = we Wor, Won, | (35) 
Wet Wore Woort 


und W(w,,, W,,) nach (33). Hierin sind w,, und w,, die beiden als bekannt anzusekende Funk- 
tionen, w die abhangige Variable, bei denen es, wie oben erlautert, hier auf die Abhangigkeit 
von t ankommt. 

Entwickelt man (35) nach der ersten Spalte, so erhalt man (34) in der Form 


w,— Aw,+ Bw=0, (36) 


worin die .,Koeffizienten“ 


ie ce a 
W (Wor We2) | Woyy, — Wong, | 
und | (37) 
pa 1 | Wor, Woo, | | 
W (wo1, Wg») Worst Woors ) 


bekannte Funktionen sind. 
Nun gilt aber auch die urspriingliche Differentialgleichung des hier behandelten 
Problems 


Wy +AAw +8 (t)w,, = 0. (38) 
Den beiden Gleichungen (36) und (38) wird Geniige getan mit 
A= 0 (39) 
und 
AAw + 0(t)w,,= Bu. (40) 


Ware A +0, so kénnte man durch Subtraktion aus (36) und (38) eine Differentialgleichung 
erster Ordnung bilden, deren Lésungen u. a. eine lineare Komposition der unabhangigen Funk- 
tionen W,, und wp. mit beliebigen Konstanten enthalten miiften, so daB z. B. das Anfangs- 
problem mit gegebener Lage und Geschwindigkeit eines Plattenpunktes lésbar ware. Das ist 
ein Widerspruch, also ist A= 0. 

Die beiden Funktionen w,, und wm, bilden also hinsichtlich ihrer Zeitabhan- 
gigkeit ein Hauptsystem fiir eine lineare homogene Differentialgleichung 
zweiter Ordnung, in welcher das Glied mit dem ersten Differentialquotienten 
nicht vorhanden ist [namlich (36) mit A= 0}. 

Hieraus folgt nach bekannten Satzen aus der Theorie der Differentialgleichungen 


W (we, We,) = konst. (41) 
und hieraus wieder, daB diese Determinante auch beim Ersatz von t durch t-+P ihren Wert 
nicht andert: 

W (wo, ? Wo,) = W (wo, > Wo,) : (42) 
Andererseits ist ; 
W (wo,, Wo,) = W [ (441 We, + 442 Woy) ’ (a5, Wo, + a9 W,)] ? 


oder wenn man die rechte Seite hiervon ausrechnet, 


Head te | y 
W (Wo, Wo,) = | 912 912! W (we, Wo,)- (43) 


| 54 G9 


1 §, hierzu: Kamke, Differentialgleichungen reeller Funktionen, Sy os UBD: 
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Aus (42) und (43) folgt die Behauptung 
=a. (44) 


g) Die fiir das Stabilitatskriterium maBbgebenden Werte bis einschlieBlich, 
Um einen Einblick in die Brauchbarkeit der Methode und die Art des Stabilitats-Kriteriums zu 
erhalten, wurde die Rechnung allgemein bis zu den Gliedern mit #4 einschlieBlich durchgefiihrt, 
indem die zu U, zugehérigen Zeitfunktionen T,, bzw. T,, nach der Darstellung (29) ermittelt, 
darin t durch t-+-P ersetzt und durch den Vergleich die Substitutionskoeffizienten bestimmt 
wurden. Das Ergebnis ist folgendes: 

Mit den Werten ,d,, und ,b,,, deren Bestimmung weiter unten dargelegt wird, sind folgende 
Ziffern definiert: 
odo 
021 Cea aaa 


Q 


b 2G 7) 
4 v= Qo adhe ees —27-2, 
Ce oul Teds Wes. bos waele 7 | 


_— 2 
= 2 (adel 5g gy = ade 93 (11 4,93), 


Wie ime: ® Uo 
Hiermit ergeben sich die zu U, gehérenden Zeitfunktionen T,, und T,, zu 
Tp, = sin Wgt + (Bo — do Yo) COS Wet — Yo Wet? sin@gt+--*, | (46) 
T o9 = COS Wo t — (Bp — Oo Yo) t SiN Wot — Po We t® cos M@et+-°-, | 


ferner die Substitutionskoeffizienten zu 
1, = (1 — Yo We PV C08 Va (Bo == 0% Yo) Psin 7, 
Ayo = (1 — Yq We P?) sin Yo +- (Bo — Oo Ye) Pos Xo » 


} 
| 
Los = Ways ( ow 
Ae (41149 — 44, = 1) , | 
oy EEN 
und schlieBlich das Stabilitatskriterium in der Form 
w beschrankt fiir t—> co wenn |cos 2 7 ho| =| (1 — a4) cos yp — ag sin y,| <1, (48) 
wobei cos 2 7 h, stets reell. 


Hierin sind noch ,d,, und ,b, zu definieren, die sich in folgender Weise ergeben: 
l. Es sei 

D i 

ah Nose (oMiixy — oftizxy) » 


wobei die Koordinatenfunktionen ,u,, und ,u,, zu bestimmen sind aus den inhomogenen Diffe- 
rentialgleichungen 


A A uy, a OoUrny 7 (@ s Wo)” in = —— 


A A ly, + Ootirexy — (© +o) oy, = + on Uoxy : 


Dann werde ,D, nach Eigenfunktionen U,, entwickelt : 


Cc 
oD, —s 2 et U,, : 
ia ——* 
Hiermit ist 
: F490 = fly, fir w=o (49) 
bestimmt. 
2. Es sei 
V 1 
OS mee (eVsexy oa oUsixy) 2 
wobei die Koordinatenfunktionen eV3, und ,v,, zu bestimmen sind aus 
AAgvs, + OnVsrny = (@e — 0)? 31 = 2(@— — @)oUyy » 


AAgv35 =+ OWsaxy — (Wp +)? og = 2 (We +O) oly, - 
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Dann sei .V; nach Eigenfunktionen U,, entwickelt: 
Vs= > eb, U, . 
“ 


Hiermit ist 
obo = ob, fir w=e@ (50) 
bestimmt. 

Damit ist das Problem mit der angegebenen Genauigkeit im wesentlichen gelést. Die ziffern- 
maBige Durchrechnung bis zum Endergebnis ist naturgema®8 umstandlich; die in den maB- 
gebenden Ziffern (45) vorkommenden Werte ,d,. und ,b, sind komplizierte Funktionen von @ 
und @, mit unendlich vielen Polen und Nullstellen, die cin besonderes Studium erfordern. 

h) Ein Beispiel in erster Annaherung. Beschrankt man sich auf Glieder mit héchstens 
#2 und legt man als ein méglichst einfaches Beispiel die quadratische Platte mit 


O= 0. | U== Py sin wt (51) 
zugrunde, so vereinfacht sich alles erheblich. Die Eigenfunktionen U, sind jetzt diejenigen 
der Platte ohne Schubbelastung am Rande, also die in Mitteilung I bereits benutzten Sinus- 
produkte. Fiir die niedrigste Eigenfunktion 

Ce iltxciny amit A,= 94, . Wp== 2 (52) 


seien die Hauptergebnisse mit einer Skizze der dazu fiihrenden Rechnungen und Uberlegungen 
kurz angegeben. 
Nach den Ausfiihrungen in Mitteilung I wird 


py 
7? oe (ux? — 1) (® —1) 
Hy? 


Hiermit sind die Differentialgleichungen fiir ,u,, und ,u,. zu lésen, die im vorigen Abschnitt 
angegeben sind und mit 0 = 0 jetzt folgendermafen lauten: 


sinuuxsin vy (u,v = 2, 4, 6,...). (53) 


Up... = C08 4 COSY = 


— (@g— @") py, +4A uyy=— 7 Ue, > 
(54) 
— (@p +)? ety, 4A Quy = airy Uo,, eo 
Mit 
i = DY %y,, Sine xsinyvy, } 
lip = Dd) %,,, Sin wx sin Vy 
by 
erhalt man die Gleichungen fiir die « und hieraus 
ia 8 py : antl 
“hus ~~ a8 [(u? £9)? — (9 — 7] (@# — 1) (*® — 1)’ | (56) 
aie 8 py | 
Lies me [(u® + 9)? — (we + w)*] (x? — 1) 0 — 1)” 
Die fir .D, maBgebende Funktion ,u,,— oly, sei angesetzt zu 
0411 — og = DS Puy sinxsin vy, (57) 
uy 
dann wird mit einfacher Zwischenrechnung 
ROPE 16 py 1 ) 
“ a —— 1) (92 Nh) eee | (58) 
ayy 
mit yy = (vw? +9)? — (w? +0”). | 
Mit (57) und (58) wird dann 
(o11 — ollya)xy = 2) Bur UY Cos fee cos VY, 2) 


Ly? 


das wieder nach den Eigenfunktionen zu entwickeln ist. 
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a 


Da w und » gerade Zahlen sind, treten jetzt in den Entwicklungen der Cosinus die ungeraden 


Vielfachen r,o auf; es wird 


(oUi¢— oU42)xy = >) Bix by ar 5 sin yx sin oy (r,¢ = 1,3, 5, ...)., (G9) 


Der Index r ist hier an Stelle des frither benutzten 9 verwendet, um eine Verwechslung mit dem 
hier anderweitig gebrauchten Index g zu vermeiden. 

Damit ist das wesentliche bereits erreicht. Nach dem vorigen Abschnitt ist, wenn hier zur 
Vermeidung von Verwechslungen ein anderer Summationsindex benutzt wird, 


1 fos} 
oD, = = (oUy1 — oW'y9) xy = & oti U, 
P= 


Der gesuchte Koeffizient ,d,, ist hierin derjenige fiir p=, die zugrundegelegte Eigenfunktion 
mit dem Index 9 ist aber in diesem Beispiel zu sin x sin y angenommen; zur Ermittlung von od1¢ 
ist also in (60) r—=o= 1 zu setzen, und man erhilt fiir das gewahlte Beispiel 


ee py 
oly. = Q Pur (uw? —1) (v2 — 1) 


oder 2 
ee: wy u 
0449 = — eine l@apeaD . __ 40% a ce) 


we Fs 


mit 
A, y= (uw? + v?)? — (wi +?) (u, v= 2, 4,6...). 
Bei der Auswertung ist zu beachten, daB a,,— a,, usw., so daB die entsprechenden Glieder 
bei « ~v doppelt auftreten. 
Damit ist der fiir das Kriterium (48) notwendige Wert ,d,, explizite gewonnen. Dieses selbst 
vereinfacht sich bei Beschrankung auf Glieder mit héchstens #? erheblich. Mit 


dl 
¢ (w) = — Wiydige 
mY 47 
X =te=2n = —, (62) 
fee 


— ic s 
und nach genauer, grundsatzlich einfacher Analyse der verschiedenen Méglichkeiten erhalt man 
fiir den vorliegenden Fall aus (48) die folgenden maBgebenden Bedingungen: 
Stabilitat liegt vor, wenn I. x? < wl(w) tg y/2 , | a 
oder, wenn II. x? < — wC(o) ctg x/2 , | ey 
und zwar in dem Sinne, dafs entweder das eine oder das andere Kriterium mafgebend ist; w ist 
hierin als positive GréBe einzufiihren. Da tg y/2 und ctg y/2 stets gleiches Vorzeichen haben 
und nur positive Werte von x? in Frage kommen, gilt das obere Kriterium, wenn ¢(w) gleiches, 
das untere, wenn es entgegengesetztes Vorzeichen wie tg ¥/2 bzw. ctg y/2 hat. Dabei lauft 
y/2 = 22/w@ von co bis 0, wahrend w von 0 bis oo lauft. 
Offenbar ist noch eine Untersuchung der neu auftretenden Funktion ¢(w) notwendig. Schreibt 
man ,d,. nach (61) aus, so erhadlt man 


Jia, 16 60 — w2 64 396 — w? 4 256 1020 — w? 
g ete 81 (60 — w)?— 16 w? | 2025 (396— w?)?— 16 w | 50 625 (1020 — w)?— 16 o® 


oder mit der ohne weiteres verstandlichen Abkiirzung 


= _ elie = Gn T5 - ? ¢ GD) = C 
ae =, b? — 16 (2) ya, , b, (64) 


= 2>— 16 


Von besonderer Wichtigkeit sind die Nullstellen von C; sie ergeben sich aus (64): 
C(w)= 0. fir dbs— l6wme— 0, 
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worin k eine der Ziffern n ist, d.h. 
C(a)= 0 fiir [(u? +9)? — i +4) — 16 w= 0, 


wobei uw, » die zu dem Reihengliede mit u—k gehérenden geraden Zahlen sind. Die leichte 
Rechnung liefert 


O, = Ww + vy? + 2 (4, ¥= 2, 4, 6, 8... .) (65) 
und damit folgende Werte fiir «,: 
6,10; 18,22; 30,34; 50,54; 70,74; 98,102; 


Es entstehen also Gruppen von je zwei Werten, die um 4 verschieden sind; die entsprechenden 
Werte dieser Zweiergruppen haben dann sukzessive je zweimal die Abstande 12, 20, 28, 36,... 
voneinander, wie dies in folgender Darstellung angedeutet ist: 


6,10; 18,22; 30,34; 50,54; 70,24; 98,102; 
12 12 20 20 cach pei hi 

Die Unendlichkeitsstellen von ¢(w) ergeben sich nicht in so einfacher Weise; es sind 
die Nullstellen der in (64) im Nenner stehenden Reihe. Man findet bei naherer Untersuchung, 
daB dieser Nenner 2'4,,, stets eine Nullstelle zwischen zwei aufeinanderfolgenden Nullstellen 
von ¢ hat und daselbst sein Vorzeichen wechselt. Die Lage dieser Stellen mu8B durch numerische 
Berechnung festgestellt werden; ihre Be- re 
deutung erscheint bis auf weiteres sekun- bp 
dar, da die vorliegende Untersuchung auf 
kleine 3 beschrankt ist. at 

In Abb. 1 dieser Mitteilung ist €(w) in 
Abhangigkeit von m aufgetragen. Fiir das 
weitere ist als Hauptsache festzuhalten, 
daB fiir 6 < w < o die Nullstellen in der 
dargelegten Weise auftreten, daB aber im 
Intervall 0 <w<6 keine Nullstellen 
von € vorhanden sind. 

Die in den Ungleichungen (63) auf- 
tretenden Funktionen tg und ctg haben 
ihre bekannten mit Vorzeichenwechsel ver- 
bundenen Null- bzw. Unendlichkeitsstellen. 
Da nach (62) 

4/2 =: as 


@ 
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so liegt bei wachsendem w die letzte 
Unendlichkeitsstelle fiir tg 7/2 und 
die letzte Nullstelle fiir ctg y/2 bei 
w= 4; es ist also 


co > tg 4/2 > 0| 
0 <ctg 7/2< a] 


Die fiir das Verhalten der Bewegung fiir standig wachsendes t bei von Null bis Unendlich 
wachsenden, jeweils konstant gehaltenen Werten von w maBgebenden Kriterien kénnen 
jetzt vollstandig diskutiert werden. 

In einer Koordinatenebene mit x? oder mit x als Ordinate und mit w als Abszisse werden 
durch die Kriterien Gebiete der Stabilitat abgegrenzt. Wahlt man zunachst fiir die Darstellung x’, 
das hier aus physikalischen Griinden nur positive Werte einschlieBlich Null haben kann, so liegen 
die Gebiete der Stabilitat in der oberen Halbebene und kénnen sich nach unten héchstens bis 
zur w-Achse erstrecken. Man stellt leicht fest, daB die Funktionen tg, ctg und ¢ in den jeweils 
in Frage kommenden Intervallen von @ monoton verlaufen; x* als Funktion von @ ist also auf 
der Grenzkurve bei diesem Grad der Annaherung stets eindeutig. Die w-Achse stellt daher 
stets die untere Grenze der Stabilitatsgebiete dar, und nach oben werden diese begrenzt durch 
die oberen Grenzkurven mit den Gleichungen 


fir 4<0m< co. (66) 


Abb. 1. Die Funktionen ¢(m) fiir das Beispiel. 
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x2 = wl(w) tg /2 im Falle I, 
v—=—wl(o)tgy/2 > 5 i. 


Dabei schlieBen sich, wie man ohne weitere Erlauterung leicht findet, die Intervalle von w 
nach Fall I von (63) mit denjenigen nach Fall II liickenlos aneinander; die w-Achse x*°=0 
ist also durchweg die untere Grenzkurve fiir die Stabilitatsgebiete, was hier aus physika- 
lischen Griinden fast selbstverstindlich ist, da bei #,—0 lediglich die freie periodische Eigen- 
schwingung der Platte ohne jede aufere Einwirkung vorliegt. 

Funktionentheoretisch sind die Nullstellen der oberen Grenzkurven, also die Stellen fiir 
Ctgy/2—0 und fiir Cctg 7/20 als Verzweigungspunkte der Kurve anzusehen, die sich aus 
simtlichen Grenzpunkten des Stabilitatsgebietes ergibt. 

Fiir die Lage dieser Punkte ist im Intervall 6 <w< die Funktion ¢ maBgebend, 
deren Nullstellen in den erwahnten, aus je zwei um 1@ = 4 voneinander entfernten Werten 
bestehenden Zweiergruppen liegen, wobei der Abstand zweier sukzessiver Gruppen bei unbe- 
schrainkt wachsendem @ ebenfalls unbeschrankt anwachst. 

Dagegen sind im Intervall 0<wm< 6 die Nullstellen von tg /2 bzw. ctg 7/2 maBgebend. 
Diese liegen offenbar bei 


(67) 


22 Tt 
fees Sikes (ite Une mas sone 
d.h. bei @= 4/m, wobei 0 < wm < 4 fir 09 Sm=Z1, wahrend das ganze Intervall 4d w< 
dem zwischen m= 1 und m= 0 liegenden Intervall von y/2 entspricht. 

Geht also@ von dem Wert 4 nach Null, so durchlauft m die ganzen Zahlen von 1 bis oo, und 
die Verzweigungsstellen liegen fiir dieses Intervall in unendlicher Anzahl bei 


O43 22 ose ls Alas 416s Ata oe: 


mit nach Null abnehmendem Abstand voneinander: Der Nullpunktist eine Haiufungsstelle 
der Verzweigungspunkte. 


Fiir die praktische Anwendung empfiehlt sich die Auftragung von x == + |x, da negative x 
ebenfalls physikalischen Sinn haben. In dieser Darstellung liegen die Gebiete der Stabilitat auf 
beiden Seiten der w-Achse. Die Grenzkurven schneiden jetzt die wm-Achse mit ,,vertikaler“ 
Tangente. Die Gebiete der Labilitat liegen ebenfalls zu beiden Seiten der m-Achse und stofen 
an diese in den Verzweigungspunkten unter Bildung von spitzen Zwickeln. Im Gebiet 0< w <4 
werden diese durch die senkrechten Geraden durch die Verzweigungspunkte und die nach links 
gekriimmten Anfange der im wesentlichen durch den Verlauf der von Null nach links ansteigen- 
den tg- und ctg-Funktionen gebildet, in dem Intervall 6 <@ < oo durch je zwei nach den 
beiden Seiten entgegengesetzt gekriimmte Kurven, deren Verlauf im wesentlichen durch den- 
jenigen von €(m) gegeben ist. 

Fiir die Anwendung ist der Kriimmungsradius der gekriimmten Grenzkurven an den 
Verzweigungsstellen von Wichtigkeit, der sich leicht berechnen 1aBt. 

Ist w, der Wert von w fiir eine Verzweigungsstelle und 6q@ eine von erster Ordnung kleine 
Abweichung von w gegeniiber @,, so daB w= wm, + 6a, so laBt sich x? leicht als linearer Aus- 
druck in der Veranderlichen dw herstellen, durch den dann in der Umgebung von «, als quadra- 
tische Funktion von x gegeben ist. Die Kriimmungsradien 9, der Kurven x(w) an den Ver- 
zweigungsstellen erhalt man dann leicht als reziproke Werte der zweiten Ableitungen von 6@ 
nach x. Die Rechnung muf fiir die beiden Intervalle 0<w< 6 und 6m < < © in verschie- 
dener Weise durchgefiihrt werden; fiir das zweitgenannte Intervall erhalt man offenbar fiir jede 
Stelle zwei Werte, die in ohne weiteres verstandlicher Weise mit Okrechts UNA Opfjinks bezeichnet 
werden. 

Man erhalt 


oe oY fiir 0< wr.< 6, (68) 


Or. 2 22 
Oklinks = 2 a, OES ae 


fir 6<a,<0, (69) 


ORF 2 Qn 
Okrechts = 2 a BO Se Sg | 
ah WO, 


wobei in (69) die oberen Vorzeichen fiir den kleineren Wert von c, fiir jede Zweiergruppe der 
Nullstellen zu nehmen sind, die unteren Vorzeichen fiir den gré®eren jw,-Wert der entsprechenden 
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Zweiergruppe. Ferner ist 


a Ey 
ESL) G2 SSN 
wobei 4 und y die beiden fiir die jeweilige Zweiergruppe von Nullstellen mafgebenden geraden 
Zahlen sind. 

Aus (68) ersieht man, daB 0, — oo fiir wm — 0; die sich nach dem Nullpunkt hin haufenden 
Stabilitatsgebiete gehen bei Annaherung an den Nullpunkt in Parallelstreifen iiber, so daB in 
der Umgebung des Nullpunktes unbeschrankt wachsende 0, erforderlich waren, um das Stabili- 
tatsgebiet zu iiberschreiten. 

Kin einfacher Uberschlag zeigt, daB die Kriimmungsradien mit unbeschrinkt wachsendem «, 
ebenfalls nach unendlich gehen. Die kleinsten Kriimmungsradien treten bei m,= 6 
und bei w,=4 auf; dies sind also die Stellen, an welchen die Gefahr anwach- 
sender Stérungen bei Vergréferung von % am griBten ist. 

In unmittelbarer Umgebung von m— ay, ist fiir die Grenzkurve 


(70) 


‘Le pre rne (71) 


worin die Betrage von 0, und dq einzusetzen sind; aus dieser Beziehung liest man ab, bis zu 


welchem x= 2 %,//~ man bei einem gewissen 6@ gehen kann, ohne da® Labilitat eintritt. 
Andererseits ist natiirlich auch 


x 
OD) == . 
20% 


(72) 


Sind also bei einer Anordnung gewisse Amplituden x, einer periodischen Belastung auf Schub 
zu befiirchten, so gibt (72) fiir jedes «, die hier vorhandene Breite dq fiir ein gefahrliches Intervall 
von j, innerhalb dessen anwachsende Stérungen zu befiirchten sind. 

Tabelle 1 gibt die einzelnen kritischen Stellen x) bei 6m = 0,1 = @,/20; in Tabelle 2 ist fiir 
dieselben Stellen dw fiir x,—10 angegeben, was ungefahr die Halfte ist von derjenigen Rand- 
belastung auf Schub, bei welcher die quadratische Platte bei den Navierschen Randbedingungen 
statisch ausbeult. 


Tabelle 1. 
AM 0667... 0:8-.| 1,0. 4.15333. 2 4 6 10 18 22 
| 
ee! 1358 | 1i27|-892" "| 657 «47 1 147 | 430 883 6622 9585 
Ck rechts Ct ON eC Sd) oo oe ie.) oo co 140 1672 50050 113 100 
0 tinks DEL O@ =0,1 | 52,12] 15,01) 13,35) 11,46| 9,13) 5,42] 9,27; 13,29 36,38 43,79 
Xo rechts beLdw@ = 0,1 | — 5,29| 18,29 | 100,5 150.4 
Tabelle 2. 
Oe | 0,667 0,8 1,0 1,333 2 4 6 10 18 22 
oF links | 
fiir x» =10 | 09,0368 | 0,0443 | 0,0561 | 0,0761 | 0.1199 | 0,3400 | 0,1162 | 0,0566 | 0,00755 | 0,00522 
Os, rechts 
firx,=10| — = laa = as — | 0,3569 | 0,0299 | 0,000999 | 0.000442 


AbschlieBend sei noch darauf hingewiesen, daB die kritischen Werte w, und diejenigen der 
Kriimmungsradien 9, an den Verzweigungsstellen als genau anzusehen sind, da sie durch korrekte 
Grenzrechnung fiir }? — 0 ermittelt sind. Beziiglich der mit der hier durchgefiihrten Annahe- 
rung ermittelten auf der w-Achse senkrechten Asymptoten der Grenzkurven ist im Vergleich 
zu den entsprechenden Kurven der gewohnlichen Differentialgleichung vielleicht die Vermutung 
gestattet, daB sie bei Beriicksichtigung geniigend vieler héherer Potenzen von J, in gekriimmte 
Kurven tibergehen. Auch sei daran erinnert, daB zu jeder Eigenfrequenz @, der unbelasteten 
Platte und der entsprechenden Eigenfunktion Uz eine andere Funktion ¢() und damit andere 
Grenzkurven gehiren. Wenn also bei der entsprechenden gewohnlichen Differentialgleichung 
die graphische Darstellung haufig als ,,Karte‘ bezeichnet wird, so kénnte man bei der hier 
behandelten partiellen Differentialgleichung von einem ,,Atlas“* solcher Karten sprechen, und 
zwar mit unendlich vielen Blattern. 


Ingenieur-Archiv 
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Abb. 2 dieser Mitteilung zeigt andeutungsweise das erste dieser Blatter, hervorgehend aus 
der hier annaherungsweise behandelten tiefsten Eigenfrequenz der quadratischen Platte. 

Durch dieses Beispiel ist, wenn auch zunachst fiir den einfachsten Fall und fiir die einfachsten 
Randbedingungen, die Anwendbarkeit der angegebenen Methoden erwiesen; die Ergebnisse 
kénnen bei einer um Null pulsierenden Schubbelastung technisch unmittelbar benutzt werden. 


2. Andere Lésungsmethoden; Ver- 
gleich der Ergebnisse. a) Ein ande- 
rer Ansatz zur Lésung. Ohne 
Ausfithrung von Einzelheiten soll in 
diesem Kapitel eine andere Me- 
thode zur Lésung der grund- 
legenden Differentialglei- 
chung entwickelt werden, die zu 
besonderen Gesichtspunkten und 
Vergleichen fihrt. 

Wesentlich fiir die Durchfiihrbar- 
keit des Potenzreihen-Verfahrens 
war die Entwicklung der in den je- 
weiligen inhomogenen Gleichungen 
auftretenden Funktionen U,, bzw. 
u,y in Fourier-Reihen nach den aus 
der ersten homogenen Gleichung 
folgenden KEigenfunktionen. Hier- 
durch wird das Verfahren gewisser- 
mafen in FluB gebracht. Dies legt den Gedanken nahe, eine derartige Entwicklung 


in dem Gliede w,y von vornherein zu versuchen. Danach wird von dem Ansatz 


J 


Vis Md YH Wy 
Yi, Gebiele anwachsender Storungen 
Abb. 2. Stabilitatskurven fiir das Beispiel. 


ee KT, (ESM (73) 
n=0 


ausgegangen, und zwar soll die gesamte Reihe der Differentialgleichung geniigen, nicht das 
einzelne Produkt K,,T,. 
Man erhalt durch Einsetzen in die Differentialgleichung 


>» T, K,+ > T, AAK, +8 > T, kK, = 0. (74) 
Ferner wird angenommen, dab Kn, nach den K, entwickelbar ist: 


ES ae wie oils 5 (75) 


=0 


womit die Summe des letzten Gliedes von (74) iibergeht in 


ASAE PR hdl ole! Pa bad (76) 


Hierin werden in iiblicher Weise unter Voraussetzung der Zulassigkeit die Indizes und die Reihen- 
folge der Summationen vertauscht; so entsteht 


R=) KS opal (77) 
n Q 
Hiermit entsteht aus (74) 
= Ty Ky ed, Pe AA Re Shoe ace) (0. ele ee Sie (78) 
n n e 


Diese Gleichung wird erfiillt, wenn fiir jedes n einzeln die Gleichung 
T, K,+T, 44K, + 0K, > an T,= 0 (79) 


oder 


(To OS oT) Ky eT eA 0 (80) 
Q 
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oder 


T’ +02 1p 
eee) | AA. 
Te ram x. a as Tame (81) 


erfillt wird, worin I, eine Konstante bedeutet. Aus (81) entstehen dann die Differentialgleichungen 
fiir die T und die U 
Tr +17, +0 > don T= 0 | 
; (82) 
AAK, —1,K,= 0. J 
Beziiglich der K, und 1, ist also eine groBe Vereinfachung eingetreten; diese sind jetzt die 
Eigenfunktionen bzw. Eigenwerte der Differentialgleichung fiir die schubfreie Platte. 
Dafiir enthalt die Differentialgleichung fiir jedes T,, simtliche T,; man hat es also nicht 


mit der einzelnen Gleichung in (82) zu tun, sondern mit einem System von unendlich vielen 
Differentialgleichungen zweiter Ordnung: 


T 7 41,Tt-8. 5) gn Tp =.0 (ae (83) 
p—v 


und 


das sich auch schreiben 148t als 
TATU oh) TLE OOS! oT, = 0: (84) 
Q 


Der Strich am Summenzeichen gibt an, daB das Glied mit p=nin der Summe nicht vorkommt. 
Ausgeschrieben lautet (84) mit 


BiamO Wen (85) 

To + (lo+ Boo) To Se ig h y + Boo Ts Ly Poy fips gee 
Bor Tot TY + (41683) Ti + Bs, Ts Ba bet 0) 

Bos Ty + Br» T,+ Ty + (g++ Bos) Te + Bye T3-++: 7 (ie 
Pos To +Bis Ty CBee TPT. 2h, 8.) Tee == 0 


h) Vergleich mit den Ergebnissen von Mitteilung I. Hierzu ist 
Ui. KONKEL. 20 0p, — konest, fur.alle 0,7 (87) 
zu setzen. Dann ist das System (86), wenn iiberhaupt, lésbar durch 
Tae Aone (88) 
g 
In bekannter Weise ergibt sich hieraus, zunadchst rein formal, mit 
rene = lire (89) 


zur Bestimmung der A, die unendliche Determinante 


I,—A B10 Boo B30 ah Vel 
Box i a Boo B31 aps amie 
D= | Boo Pip LA igg ewt*s) aco * 0.5 (90) 


Bos Pas Pos LA ‘ 


4, ist eine der 9 Wurzeln dieser Gleichung fiir A. 

; Derartige unendliche Determinanten sind bekannt, ebenso die Bedingungen, unter welchen 
sie konvergieren, sowie zweckmaBige Methoden fiir die Auswertung. Jedenfalls gibt es 
einen zulassigen Bereich der maBgebenden Gréfen, fiir welchen D konvergiert 
und Lésungen hat, fiir dén also das ganze Verfahren legitimiert ist. : 

In diesem Bereich erhalt man unter Beschrankung auf einige Hauptpunkte folgenden Ver- 
gleich fiir die Ergebnisse der beiden Methoden: 
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1. In Mitteilung I ergab sich bei @ = 0) = konst. eine allgemeine Lésung in der Form 
Tien Wit We (piel ee eon ees 
19 
hierin sind die U, die Eigenfunktionen der mit J) belasteten Platte (91) a 
TA, 610 Yap t -- B, cos Was b., 
und A, der zu U, gehérige Higenwert. 


In der hier durchgefiithrten Untersuchung ergab sich dagegen 


w= SUB ES (92) 
hierin sind die K, die Eigenfunktionen der schubfreien Platte; dabei ist B. 


T,, = > (Ag sin Vagt + By cos Vat), 


Q 
und die A, ergeben sich aus D= 0. 


Bei A sind die U, angegeben als unendliche Reihen in den K,. Multipliziert 
man also (91) aus und faBt jeweils die Glieder mit gleichem K, zusammen, so erhalt man die 
Lésung in der Form B. Da hier nur eindeutige Lisungen in Frage kommen, so ergibt sich: 

Die aus (90) folgenden 4, sind die A, in A, und jedes der beiden Ergebnisse 
ist eine Umordnung des anderen. 

2. Physikalisch ergeben sich aus A und B zwei verschiedene Anschauungsweisen fiir 
die Bewegung der Platte: 

Nach A wird die Gestalt der Platte in jedem Augenblick aus den durch die U, gegebenen 
Grundformen aufgebaut, also denjenigen der mit J) 40 belasteten Platte. Fir jede dieser 
Grundformen ergibt sich die zugehérige Anderung mit der Zeit durch eine einfache harmo- 
nische Funktion mit der Kreisfrequenz |//,. Die verschiedenen A, stehen zueinander in keinem 
rationalen Verhaltnis; die Gesamtbewegung ist nichtperiodisch. 

Nach B wird die Gestalt der Platte in jedem Augenblick aus den durch die K, gegebenen 
Grundformen der schubfreien Platte aufgebaut. Fiir jede dieser Grundformen ergibt sich 
die zugehérige Anderung mit der Zeit durch ein Spektrum harmonischer Funktionen von 
unendlicher Komponentenzahl, wobei fiir jede Grundform der schubfreien Platte samt- 
liche 2, in dem Spektrum auftreten. 

3. Nach B ergeben sich die A als Wurzeln der Gleichung D0, worin D eine unendliche 
Determinante ist. 

Nach A ergibt sich jedes 4, einzeln aus einem expliziten Ausdruck fiir dieses A, allein, der 
ohne Riicksicht auf die anderen unmittelbar auszuwerten ist. 


A gibt also eine Methode zur expliziten Berechnung der einzelnen Wurzeln 
der Determinantengleichung (90), jede als Potenzreihe in %, das als gemeinsamer Faktor 
in den Koeffizienten 6 der Determinante enthalten ist. 

Auf weitere Einzelheiten sei hier nicht eingegangen. 


c) Vergleich mit Abschnitt 1 dieser Mitteilung. Ist # nicht konstant, sondern, wie in 
Abschnitt 1 dieser Mitteilung vorausgesetzt, d= O + 0) sin wt, so lautet (84) 


TY 11, T, + (O18, sin ait) Sao, To= 0 (ions eee (93) 
o=0 


Man kénnte diesen Gleichungssatz als ein System von unendlich vielen Mathieuschen 
Gleichungen bezeichnen. 

Dieses soll hier nicht direkt bearbeitet werden. Doch ist es indirekt, im Zusammenhang mit 
der Plattentheorie, in Abschnitt 1 der vorliegenden Mitteilung gelést. Wenn man in der dort 
autgestellten Lésung die U, in den K, ausdriickt, dann die Glieder mit jeweils gleichem K, 
zusammenfaht und dieses ausklammert, so stellen die Klammerinhalte Lésungen von (93) dar, 
entwickelt nach Potenzen des Parameters %o. 


(Eingegangen am 7. Juni 1950.) 


Anschrift des Verfassers: Professor W. Kucharski, Berlin-Charlottenburg, Hardenbergstr. 34. 
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». $13, 


5. 316, 


S. 316, 


SLi 


S. 318, 


5. 319, 


5. 322, 


vor der Formel in der Mitte der Seite mu® das Wort ,,Aus‘ hinzugefiigt werden. 


mu heiBen: 
Aus M=v,(0,F, +o02F.), usw. 


der letzte Teil der Gleichung (24) im Nenner muB heifen: 
1 1 
(Pa) = statt (2\* 
P De 
am Anfang der Gleichung (24a) mu es heiBen: 
OoFeve statt 0) = yx? 
2. Zeile von unten: 


d) Sonderfall (F,— Fy,)=0. statt (F,—F,)=0. 


das zweite Wurzelzeichen ist zu lang. Die Formel sieht so aus: 


ye| GF le) 


die Formel in der Mitte der Seite mu heiBen: 


ME — = Peacatatk ME- ip 


12. Zeile von unten: 


,,Beispielstrahlrohres“‘ statt ,,Beispielstahlrohres‘‘ 
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M Kinfuhrung in den Warme- und Stoffaustausch 


Inhaltstibersicht: 
I. Die Grundbegriffe des Wirmeaustausches Il, DiciiWarnielsitang. d é 
i 4 ] oo TK eleitung. — III. Der Wirmetibergang: A. Grund- 
Prete der Stromungslehre. B. Erzwungene Konvektion in lauiteret Strémung. C, Eeeohinee moneetion. oi 
urbulenter Strémung. D. Freie Konvektion. E. Kondensation und Verdampfung. — 1V. Die Wiarmestrahlung. 


A. Die [ 
ee auer eee B. Der Step blungsaustanseh. — V. Der Stoffaustausch. — Anhang. — Namen- und Sach- 


Einfuhrung in die Technische Thermodynamik 
undindie Grundlagen der chemischen Thermodynamik 
Von 


Dr.-Ing. Ernst Schmidt ; ) 


o. Professor und Direktor des Instituts fiir Wirmctechnik 
an der Technischen Hochsghule Braunschweig 


Vierte x itberarbeitete und erweiterte Auflage. Mit .244 Abbilduggen und 69 Tabellen 
im Text sowie 3 Dampftafeln als Anlage. XVI, 520 Seiten. 1950. Ganzleinen DMark 30.— 


y Inhaltstibersicht: 

I. Temperatur und Warmemenge. — II. Erster Hauptsatz der Warmelehre. — III. Der thermodynamische Zustand 
eines K6rpers. — IV. Das yollkommene Gas. — V. Kreisprozesse. — VI. Der zweite Hauptsatz der Warmelehre. — 
Vil. Anwendung der Gasgesetze und der beiden Haupisitze auf Gasmaschinen. —. VIII. Die Eigenschaften der 
Dampfe. — IX. Das Erstarren und der feste Zustand.. — X. Anwendungen auf die Dampfmaschine. — XI. Zustands- . 
gleichungen von Dampfen. — XII. Die Verbrennungserscheinungen.. — XIII. Strémende Bewegung von Gasen . 
und Dampfen. — XIV. Strémungsmaschinen. — XV. Thermodynamik des Raketenantriebes: — XVI. Thermo- 

_dynamischer Luftstrahlantrieb. — XVII. Die Grundbegriffe der Warmetibertragung. —- K VIII. Die Warme- 
ubertragung durch Strahlung. — XIX. Dampf-Gas-Gemische. — XX. Die Anwendung des I, und II. Hauptsatzes 
der Thermodynamik auf chemische Vorgange. —- XXI. Das Nernst’sche Wirmetheorem oder der dritte Hauptsatz 
der Warmelebre. — Anhang: Dampftabellen und Tafeln. — Namen- und Sachverzeichnis, 


- Die Grundlagen der angewandten 


2ae - Thermodynamik 
Von 


Dr.-Ing. habil. Kurt Nesselmann ie 
Wiesbaden 


Mit 311 Abbildungen und 5 Diagrammen im Text 
XI, 320 Seiten, 1950. Ganzleinen DMark 18.— 


y ‘ Inhaltstibersicht: . 

f I. Grundbegriffe. — II, Die vollkommenen Gase. — III. Die Dimpfe. — IV. Die strémende Bewegung der Gase 
und Dampfe.~— V. Die unvollkommenen Gase. — Vi. Thermodynemik der Gemische. — VII. Thermodynamik 
der chemischen Reaktionen. — VIII. Warmeaustausch. — 1X. Stoffaustausch. — Namen- und Sachverzeichnis. 


Leitfaden der Technischen Wirmelehre 


nebst Anwendungsbeispielen 
Von * 


Dr.Ing. habil. Hugo Richter 


: Gummersbach 
Mit 384 Abbildungen, 1 Diagramm und’ 104 Zahlentafeln 
im Text und Anhang. XII, 617 Seiten.. 1950. Ganzleinen DM 34,50 


Inhaltsttbersicht: 
A. Allgemeine Grundlagen der technischen Warmelehre, I. Allgemeine physikalische Grundlagen. — II. Feste und 
fliissige,Kérper. — III. Vollkommene Gase. — IV. Der 1. Hauptsatz der mechanischen ‘Wirmetheorie. — V. Wirk- 
liche Gase. — VI. Zustandsinderungen von Gasen. — VII. Kreisprozesse von vollkommenen Gasen. —- VIII. Der 
2. Hauptsatz der mechanischen Warmetheorie. — IX. Entropie der vollkommenen Gase. — X. Umkehrbare Kveis- 
prozesse der Gase im Warmediagramm. — XI. Entropie anderer als gasférmiger Kérper. — XII. Warmeinhalt 
der Gase bei konstantem Druck. — XIII. Nichtumkehrbare Vorgiinge. — XIV.. Dimpfe. — XV. Mischungen, — 
XVI. Strémende Bewegung von Gasen, Dimpfen und Flissigkeiten. — XVII. Warmeibertragung. — XVIII. Ver- 
brennung. — B. Die wichtigsten technischen Anwendungen. Kraftprozesse. I. Arbeitsweise der Dampfkraft- 
maschinen. — II. Arbeitsweise der Verbrennungskraftmaschinen. — III. Vergleichende Betrachtungen zu den 
P Warmekraftmaschinen. — Arbeitsprozesse. IV. Gasverdichter. — V. Kaltemaschinen. — VI. Warme- 
‘ pumpen. — VII. Gasverfliissigungsanlagen. — VIII. Anlagen zur Trocknung mit Luft. — IX, Riickkihlanlagen. — 
a. MeBverfahren. X. Messung der Luftfeuchtigkeit. — XI. DurchfluBmessung mit Drosselgerditen. — Ver- 
gasungvonfestenBrennstoffen. XII. Gasgeneratorprozesse. — Sachverzeichnis. 
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| Dr.-Ing, habil. Ernst Eckert 
Mit 125 Abbildungen. VII, 203 Seiten. 1949, DMark.21.—, Ganzleinen DMark 24,— | 
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Technische Physik | 


in Einzeldarstellungen: 
Herausgegeben yon W. Meissner 
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H ochstromkohlebogen. Physik und Technik einer Hoch- 
temperatur-Bogenentladung. Von Professor Dr. Wolfgang Finkelnburg. Mit 


132 Abbildungen, VIII, 221 Seiten, 1948. DMark 22.50 


Grundlagen der Hochstfrequenztechnik. 


Von Dr.-Ing. F. W. Gundlach, Professor an der Technischen Hochschule Darmstadt. 


Direktor des Instituts fiir Fernmeldetechnische Gerate und Anlagen. Mit 189 Ab- 


‘bildungen, VIII, 499 Seiten, 1950. - DMark 48.— 


Das voriltgevde Buch umfaBt alle wesentlichen Teilgebiete der Hochstfrcquenztechnik wie Réhbren, 
Leitungen, ife)lleitungen, Resonatoren, Schaltungselemente und Antennen. Der Titel ,,Grundlagen*™ 
ist so aufzufassen, daB® die fiir die Héchstfrequenztechnik grundlegenden Zusammenhange in wissen- 
sehaftlicher Strenge dargestellt sind. Es handelt sich nicht um’ eine erste Einfihrung in das Gebiet. 
Die Darstellung geht von den physikalischen Grundlagen (Elektronentheorie und Maxwellsche Feld- 
theorie aus und baut darauf die Erscheinyngen der Elektronenstrémung im Laufzeitgebiet und die 
Wellenausbreitung in Doppelleitern, Hohlleitern, Resonatoren und Antennen in einheitlicher Be- 
trachtungsweise auf, Hin ausfihrliches Literaturverzeichnis am SchluB® des Buches gibt die erforder- 
lichen’ Quellénnachweise und weist dariiber hinaus den Leser auf ausfithrlichere Arbeiten und auf Ge- 
biete hin, die innerhalb des gespannten Rahmens nicht gebracht werden konnten. | 


Warmeaustausch im Gegenstrom, Gleichstrom 


und Kreuzstrom. Von Dr.-Ing. Helmuth Hausen, 0. Professor an der 
Technischen Hochschule Hannover, Mit 230 Textabbildungen. XII, 464 Seiten. 
1950. See ted DMark 69,— 


Im vorliegenden Buch werden die, Gesetze der Warmeiibertragung zwischen Stoffen, die sich im Gleich- 
strom, Gegenstrom oder Kreuzstrom bewegen, yom streng physikalischen Standpunkt aus so behandelt, 
daB sie auch dem Praktiker mit hinreichendmathematisch-physikalischen Kenntnissen verstaéndlich 
erscheinen, Den Ausgangspunkt bilden die Ergebnisse der Forschung tiber den Warmeiibergang und 
Druckabfall in Rohrleitungen und Kanilen. Die Hauptaufgabe des Buches besteht darin, alle wichtigeren 
Theorien, die bisher iber Warmeaustauscher entwickelt worden sind, zusammenfassend und einheitlich 
darzustellen. Am Schlusse jeder Erérterung soll gezeigt werden, wie man die selbst aus verwickelten 
Theorien sich ergebenden Berechnungsverfahren fast immer verhaltnismaBig einfach und rasch an- 
wenden kann. \ ¢ ‘ 4 , 2 ; 


Kleinste Drucke — ihre Messung und Erzeugung. Von Dr. Rudolf 
Jaeckel, apl. Professor der Physik an der Universitat Bonn und Leiter der Hoch- 


vakuumabteilung der Firma E. Leybold’s Nachf., Kéln-Bayental,, unter Mit- 
arbeit von Dr, Helmut Sch warzund Dr, Elisabeth Schiller. Mit 301 Text- 
abbildungen. X, 302 Seiten. 1950. sec DMark 39.60 


Im Hinblick auf die in stiirmischer Entwicklung befindliche Anwendung der Hockivabunmteshatt in 


Wissenschaft und Industrie nicht nur in Deutschland, sondern besonders auch in den USA diirfte die- 


Darstellung ihrer physikalisch-technischen Grundlagen durch den Verfasser nicht nur fiir Physiker, 
sondern auch fir den auf den Nachbargebieten Arbeitenden (Chemiker, Pharmazeuten, Mediziner, In- 
genieur usw.) gerade jetzt von#akutem Interesse sein. Der Verfasser hat als Leiter der Entwicklungs- 


abteilung einer auf dem Gebiet der Hochvakuumtechnik fihrenden Spezialfirma und als letzter Mit- _ 


arbeiter Gaedes besonders enge Berthrung mit allen im Vordergrunde des Interesses stehenden 
Problemen. Bay it e 


Die elektromagnetische Schirmung in der Fern- 
melde- und Hochfrequenztechnik. von pr. phil. Heinrich 
Kaden, Ober-Ing. der Siemens & Halske A.-G. Mit 145 Textabbildungen., VIII, 
274 Seiten. 1950. DMark 38,— 


Ein aus der Forderung nach Wirtschaftlichkeit geborener Grundsatz in der Technik ist der, eine ver- 
langte Wirkung mit méglichst geringem Aufwand zuerzielen. Zur Verwirklichung dieses Prinzips muS 
sich die Technik der Methoden und Ergebnisse der exakten Naturwissenschaften und der Mathematik 
bedienen. Je mehr die wissenschaftlich arbeitenden Ingexieure hiervon beherrschen, um so wirtschaft- 
licher und exakter werden sie ihre technischen Aufgaben lésen kénnen. In diesem Sinne soll ihnen 
dieses Buch Unterlagen auf dem Gebiete der elektrognagnetischen Schirmung in die Hand geben, Das 
Buch enthalt viele Ergebnisse, die bisher noch nights Goabtiesent worden sind und daher fiir die breite 
Fachwelt als neu bezeichnet werden kénnen. : é { 
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